


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1981 


Time domain radar laboratory operating 
system development and transient EM analysis. 


Sorrentino, Ludwig A. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/20696 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








DUDLEY KNox LIBRARY 
NAVAL POSTGRADUATE 


SCHOOL 
MONTEREY, CALIF, 936 


40 





MPG - i7Y 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





imeSIsS 


TIME DOMAIN RADAR LABORATORY OPERATING 
Seo Pee DE VELOEMENT AND TRANSIENT 
EM ANALYSIS 


by 


Pua? Pe SOrrentino 


september 1981 





Thesis Advisor: M. A. Morgan 


Aporoved fOr public release; distribution unlimited 


T202794 











SECURITY CLASSIFICATION OF THIS PAGE (When Dette Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO. 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


S$. TYPE OF REPORT & PERIOO COVEREO 
Master's Thesis; 








TITLE (and Suatitie) 


TIME DOMAIN RADAR LABORATORY OPERATING SYSTEM 





4. 
















DEVELOPMENT AND TRANSIENT EM ANALYSIS ekaSiilole ae 
- AUTHOR a) 8. CONTRACT OR GRANT NUMBER(8) 





Ludwig A. Sorrentino 


A ADORESS 10. PROGRAM ELEMENT, PQOJECT TASK 
. PERFORMING ORGANIZATION NAME ANO AO Rae cnee eGo ENtae ONES: 


Naval Postgraduate School 
Monterey, California 93940 


CONTROLLING OFFICE NAME ANO AOORESS 12. REPORT OATE 


198 
Naval Postgraduate School ~peptemer Set 


Monterey, California 93940 298 


. MONITORING AGENCY NAME & ADORESS/(If different trom Cantraiiing Ollice) 18. SECURITY CLASS. (of tAte report) 


Naval Postgraduate School 


Monterey, California 93940 Unclassified 


Se. OECLASS!I FIC ATION/ DOWNGRADING 
SCHEOULE 


16. OLSTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited 


. DISTRIBUTION STATEMENT (al the adetrect entered in Block 20, if different from Report) 








. SUPPLEMENTARY NOTES 






. KEY WOROS (Continue on reverse side if neceseary and identify by block 2umber) 


Transient imaging, ramp response synthesis, transient electromagnetics, 
electromagnetic measurements, electromagnetic scattering, radar target 
identification. 







. ABSTRACT (Continue an reveree side if neceesary and identity by blieck mamber) 


An operator-microprocessor interactive Operating System has been developed 
for the Time Domain Radar Laboratory (TDRL). The Operating System performs 
Signal acquisition and averaging, real time and frequency domain computa- 
tions and provides outputs in easily evaluated graphic displays. Target 
Classification is made by analysis of either impulse, step or ramp responses. 
A nolise-reduction signal optimization technique is implemented. Numerous 
measurements of known and unknown targets are made using various antennas 









DO ' hag 1473 EOITION OF ! NOV 68 18 OCBSOLETE 


JAmM 73 
Je - —— ee 
S/N 0107-014- 6601 | SECURITY CLASSIFICATION OF THIS PAGE (When Dere Entered) 








ER NE ne 
Fecumsty CL AGEIFIC ATION OF Twit PAGE Mere Enz >), 






and results are compared to theory. Targets are classified. Antenna 
parameters are established. Algorithms using the specialized features 
of the host Graphic System are tailored to the requirements of the TDRL 
for a detailed graphic presentation of processed data. 


DD orm. 1473 
o/N ‘01922014-6601 


ee nee 
SECUAITY CLASSIFICATION OF Tris PacerBren Cate Entered) 








Approved Zor public release; distribution unlimited 


Time Domain Radar Laboratory Operating System Development 


and Transieat iM Analysis 


Dy 


Ludwig A. Sorreatino 
Lieutenant Commander, United States Navy 
B.S.O.E., United States Naval Academy, 1971 


Submitted in partial fulfillment of the 


requirements for the degree of 
Dae aeor SCLENCE IN SLECTRECAL ENGINEERING 


NAVAL POSTGRADUATE SCHOOL 
September 1981 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIF, 93940 


ABSTRACT 
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Algorithms using the specialized features of the host 


Graphic System are tailored to the requirements of the TDRL 


for a detailed graphic presentation of processed data. 
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Le. INTRODUCTION 


A. OVERVIEW 

Seaetly stated, direct time-domain technigues for 
transient cadiation and scattering are broadly concerned 
With the response of some object of interest, either 
equipment or a target, +t95 non-monochromatic electromangetic 
sources. The analysis is done in the tine-domain, with time 
treated as an explicit independant variable in_ most methods. 

In the past, frequency domain analysis has been the 
primary means of analyzing non-monochromatic EM phenomena. 
Time-domain solutions in closed-form wera almost impossible 
to obtain, except in certain very Simple cases. The amount 
of computational work reguired was correspondingly 
Significant. No practical aides were available to assist in 
performing the operations prior to the development of fast 
computer systems. 

ApoOerege Limitation was the requirement that the 
time-domain response be as nearly an impulse as possible and 
thus possess a broad band of frequencies. fardware was 


reguired to produce very short pulse waveforms. 


16 





The situation changed radically with the widescale 
Mesoates2.on Of computers capable G2" Jo beonland S rapid 
computations involving large numbers Of pe terattons |) ~2n 
relatively short periods. pi coupled with the rapid 
maturing of short-pulse hardware technology, soon made 
possible the application of time-domain techniques to a wide 
variety of eletromagnetic problems. 
Miller and Landt [Ref. 1} list some advantages of 
direct tine-domain solutions over frequency domain 
treatments of transient problems. They are: 
1. Greater solution efficiency for many types of EM 
problems. 

Zeemote Convenient handling of non-linearities. 

3. Improved physical insight and interoretablility of EM 
phenomenon. 

4. Availability of wide-band information from a single 
Calculation. 

5. Opportunity to isolate interactions, using tine range 
gating (e.g., pulse reflections fron distant objects 
Nomewon taterest, etc.,). 

The major disadvantage i the increased complexity of 
the computer implementation of a time-domain program as 


opposed to one involving freguency domain techniques. There 
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is also a corresponding GEELLCULty ‘in the software 
development and its Practreal use. Famally, the 
computational time may be significantly greater than that 
required of other techniques. 

The applications for time-domain techniques are diverse. 
Transient analysis can be used to directly evaluate the 
transient charact2ristics oie new antenna designs, 
determining such parameters as driving point currents, 
radiated far fields and near fields, and assisting in 
determining antenna gain. These methods may also be used in 
determining the degree of coupling due to electromagnetic 
pulses (EMP) generated because of lightening during a 
nuclear detonation. Finally, time-domain technigues may be 
directly applied in the solution of the inverse scattering 
problem involving the receiving and scattering properties of 
a large category of target types. 

As time-domain transient techniques are further refined, 
it can be expected that a greater number of uses will be 


identified. 


Be. REVIEW OF TRANSIENT EM ANALYSIS TECHNIQUES 
There are three broad categories of techniques 


identifiable in the literature involving transient analysis. 
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They are direct (time-donain) techniques, transform 
(frequency domain) techniques, and hybrid 
wecnmaques. [Ref. 2] Miesee 6rSpOrteis only winterested™ in 
direct techniques and will not review the others. 
lo memvyesGel Optics Inverse Scattering 

The earliest technigues developed involved physical 
optics inverse scattering. Kennaugh and Cosgriff [Ref. 3] 
in 1958 proposed that the impulse response of a scattering, 
conducting body is simply the second darivative of the 
projected area function of the body if physical optics 
currents are postulated on the surface. Ifa target is 
illuminated with an incident ramp plane wave, then the 
projected area function of the target is directly 
proportional to the far-scattered ramp response. This is 
only an approximate relationship, however, since it is based 
on the assumption of physical optics currents. 

A more direct and precise ramp response technique 
WaS suggested by Kennaugh and Moffat (Ref. 4] in 1965. 
Mhey predicted Satisfactory apvroximate ramp response 
Signatures could b2 obtained with NMacrower bandwidth 
interrogating signals than that required for the impulse 


response. 
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The above technigues were generally applicable for 
targets which were rotationally symmetric and for whom the 
direction of the electromagnetic field incidence was axial. 
In 1976, Young { Ref. 5] demonstrated that a volume estimate 
can be obtained from the ramp response waveform. He further 
extended the ramp response technique to arbitrarily shaped 
targets with the incident field being other than axial. 

eee ca cel Tnverse Scattering Method 

Another method which improves upon the physical 
optics inverse scattering theory is the "exact" inverse 
scattering method first suggested by Bennett et. oy 
meeremo}  2n 1974, and further refined in 1977. [kRef. 7] 
The inverse scattering problem is solved using an inversion 
of a space-time integral. The method takes into account the 
backscatter response of the target, the direct effect of the 
incident electromagnetic field, and the correction currents 
flowing on the target surface. Theoretically, by completely 
describing the interactions of the two components of the 
surface currents, a close approximation of the target shape 
is possible. The actual implementation is limited to 
Meeriteccriy"” conducting scatterers. * is this method which 
has been selected for use in the Naval Postgraduate School 


Time Domain Radar Laboratory (TDRL). 
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3. Other Methods 


Other methods that have been developed, although to 
a less satisfactory degree of accuracy than the inverse 
scattering methods, include radar target identification 
methods. Only a £inite number of solutions (targets shapes) 
are allowed as opposed to the nearly infinite number 
possible using inverse scattering. The models of targets 
are stored for parameter comparison with actual ‘ceturns. 
Certain parametric aspects of a wideband illuminating 
signal, such as natural resonance frequencies and 
polarization of the return, are used. Other similar methods 
use amplitude and phase differences or special properties of 
Rayleigh region scattering. All these methods are linited 
in practical application. 

Sty JL another important and powerful method 
considers multifreguency target illumination using existing 
radars to simultaneously illuminate targets with VHF, UHF, L 
and S band radar frequencies with circular polarization. 
{[Ref,. 8} Circular polarization eliminates any polarization 
angle dependance the signal might have. This approach 
aliows for the inclusion of ail radar system parameters 
(such as antenna gain, noise level, range, etc.) in the 


maent. recat lon algorithna. The technique [ue emer 


Zul 








demonstrates that the target shane information may be 
obtained with a high degree of accuracy from radar systems 
Genrectlys tailored to the probiem of general target 


moeentification. 


C. REVIEW OF TDRL PROGRESS AT NPS 

Three requirements were initially identified. Best ere 
method had to be chosen from the techniques listed above. 
Secondly, the construction of a physical laboratory to 
perform the required measurements in support of the chosen 
technigue had to be designed and built. Dalida an 
alogorithm had to be developed employing the selected 
technigue for laboratory use. 

AS stated above, the "exact" inverse scattering method 
suggested by Bennett was the scheme chosen. 10 SUppeEC 
this, work began in January 1980 where evaluation of various 
designs to pe used in constructing a Time Domain Radar 
Laboratory at the Naval Postgraduate School was carried out 
by Capt. C. W. Hammond, USMC, as a thesis project. work was 
essentially completed with the erection of the laboratory 
and interfacing of hardware in September, 1980. [Ref. 9] 

Inmem@opore of the TDRL, LCDR Meir Morag, Isreali Navy, 


began, as a thesis project, the development of an “exact! 


IATA 











inverse scattering algorithm in September, 1980. The 
algorithm was completed in March 1981 [Ref. 10] However, 
the program was developed in FORTRAN, a Language not used by 
the TDRL computer systen. A translation toweszasic and 
adaptation of the cod2 was necessary «95 reduce memory 
storage requiraments and to 2ansure computational times were 
of reasonable lengths. 

The final step in the establishment of the TDRL was the 
validation of the laboratory in confirming the work of 
Hammond and Morag. This is one of the principal aims of 


@hvs project. 


D. THESIS OBJECTIVES 

The objectives of this th2sis are to: 

1. Develop a working software program which will acquire 
and prepare transient response data ‘from various 
targets for uipue co the time-domain integral 
equation. 

2. To show that through signal processing, the effects 
of the receiving antenna can be neglected, and the 
response to arbitrary specified exciting waveforms 
can be constructed, Crem SyMitnesized™, wusing FFT 


technigues. 
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Bemercligate the TDRL. 

The main effort has been directed toward the development 
of a user oriented computer progran which will allow the 
determination of the second objective and prepare the 
transient domain daca, for input ints the numerical 
time-domain integral equation. The shape of the object 
being observed is displayed both before and after correction 
factors are applied. The shape of the object was developed 
from the ramp response of the impulse function applicable to 
ee target. The program was also made compatibie with the 
input of data to an alternate natural resonance frequency 
target identification technigue, using Prony's method, which 
is being developed by LT Demetrius Papaspiridakos, Hellinic 
Navy, as a thesis project. 

This report first describes the physical structure of 
the TDRL. It then derives the relevant time-domain integral 
equations for both the raiiation and inverse scattering 
solution for simple axisymmetrical, perfectly conducting, 
closed surfaces in Thapter [II. Next, it discusses and 
develops the numerical solutions to the integral equations. 
In Chapter V, the implementation of the complete computer 
code is presented. And finally, measurements that are 
relevant to the validation and confirmation of previously 


developed concepts are reported. 
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G@aes) ©Resis is the third to be produced for this major 
BEQjvect at the Naval Postgraduate School. It marks a point 
of departure in that essentially all the preliminary work 
Which will allow for the effective use of the TDRL has been 
performed with validation. The laboratory is, fer the most 
part, ready for applied research involving transient 


time-domain anaiysis. 
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ile  Paveseae DE PTION OF TODRL 


A. INTRODUCTION 

As stated in Chapter I., the primary thrust of this work 
is to provide a practical working alogoritha to serve as the 
principal link between the acgyuisition of transient data for 
targets situated on the TDRL scattering range, and the use 
of that data in the inverse scattering algorithn. 

The work performed on this project is a direct extension 
of the preliminary work performed by Hammond in the 
construction of the Naval Postgraduate School Time Domain 
Radar Laboratory (TDRL), and by Morag in the development of 
the numerical integral equation algorithm for the solution 
of the inverse scattering problem involving axisymmetric 
targets. 


The results obtained by Hammond and Morag are 


2ssentially complementary. However, the independant 
approaches each took in their work, particularly in the 
selection of unrelated processing machines for data 


processing, made necessary certain practical modifications 
to the TDRL. It is germane to review the basic work of 


Hammond and discuss those mojifications aade2. Such a review 
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Weerteasso dilow for the setting of a frame of referance to 
aid in the explanation of the linking algorithm developed 


and resuits obtained in this project. 


Be PHYSICAL DESCRIPTION OF TDRL 
The physical set-up of the ITDRL is shown in Figure 2.1. 
moecanheepe seen, the TDRL can be broken functionally into 
three groups: 
1. The Impulse Generator Source Group 
2. The Imaging Plane Group 
3. The Sampling Receiver and Signal Processing Group 
Basically, the TDRL is an open system designed to determine 
the transfer function of atarget from its backscattered 
radiation. Scattering measurements are bi-static. A 
description of each functional group follows. 
1. Source Group 
PnesmeegEOouomeansists §of a single iten, 2 Video 
Communications UHF Impulse Generator, Model 1000. This 
generator is the source of the transmitted pulse and also 
provides the triggering signal which synchronizes the Signal 
Processing Sroup to the transmitted signal. 
The Impulse Senerator is specified to provide a 


fixed amplitude 45- volt pulse with a maxinum frise tine of 
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Pigure 2.2. Block Diagram of YDRL 
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Pigure 2.3. Source Pulse from Generator 


370 picoseconds. Figure 2.3 is atypical waveform as 
measured on a Tektronix S-6 Sampler Head by the 
manufacturer. 

The output of the Impulse Generator as measured at 
the TDRL as measured by the Operating System is given in 
Peagurie 224, and Figure 2.5 as photographed from the DPO 
directly. To set full waveform presentation on the measuring 
oscilloscope, 30dB attenuation was applied to the generator 
Sut put. Taking the attenuation into consideration, the 
measured peak value of voltage (1.4604 volts) gives the 
actual output of the generator to be 46.2 volts. This is 


about 3% better than rated. The rise time of the output is 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
1 6-26-81 N/A HO 777 GENERATOR WAVEFORM 


MAXIMUM PEAK VALUE (VOLTS) 1.1284 
MINIMUM PEAK VALUE (VOLTS) -8.3320 
RMS VALUE (VOLTS) 0.4214 
MEAN VALUE (VOLTS) 0.00898 
eee or Seems os cere ae se leet VOLVER 


DIRECT WAVEFORN 
1.80 


8.58 


E-3 SEC 
0.58 1.08 1.58 2.88 


Figure 2.4. Source Waveform Measured at TDRL 


about 250 picoseconds, also better than specified by the 
Manufacturer, by about 48%. 

Some notable features of the generator used in TDRL 
measurements are a variable trigger repetition rate from 10 
GzgeeecOo 1 MHZ and an adjustable delay from 50 to 150 
NManoseconds. Ail outputs aré matched into a 50-ohm load. 
The generator allows rapid acguisition of highly coherent, 
repetitive waveforms. This 1S an important feature when 
attempting to maximiz® signal-to-noise ratio through signal 


averaging. 
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Figure 2.5. DPO Waveform--UHP Impulse Generator 


The Video Communications Impulse Generator replaces 
a Tektronix Type 109 impulse generator used by Hammond. The 
completely solid-state electronics of the Mark 1000 has 
proven more reliable in producing a relatively noiseless, 
pulse~to-pulse coherent signal than did the reed-switched 
Type 109. The Type 109 had a fixed pulse repetition rate at 
about 720 pulses per second. This is a very limiting factor 
when attempting +o noise average a signal. 

The impulse source is connected to the transmitting 
antenna via a 3.2 meter RG-8A/U coaxial cable. This cable 
serves both as a conductor of the source puise anda short 
delay element to the transmitted signal. 
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26 Imaging Plane Group 


The Imaging Plane Group consists of four elements 
the transmitting antenna, the receiving antenna, he 
target, and the imaging plane. 

a. Transmitting Antenna 

The transmitting antenna is a 6.4=-meter endfed 
nonopvole. Zt extends vertically from the ground plane and 
produces a vertically polarized EM field. 

b. Receiving Antenna -- Transverse Zlectromagnetic 

Horn (TEM) 

Hammond explored a variety of receiving antenna 
types +o be used in the IDRL. The most useful is the TEM 
horn antenna in that it has the following advantages over 
the other antennas considered: 

iepetter directivity. 

2. Less Attenuation. 

3. Mobility on the image plane. 

%. A readily variable characteristic impedance. 
5. Relatively frequency independant over a wideband. 

The TEM Horn used in the TDRL is physically an 

open sided pyrimidal structure when deployed. It is a 
single piece of flat copper, about tmm thick and trianguiar 


in shape. Figure 2.5 gives the relevant dimensions. The 
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f-——————_— 432 mm ——_——_—_——| 


Pigure 2.6. TEM Horn Dimensions 


hole at the apex of the antenna allows for connecting to the 
receiving cable. 

The receiving horn antenna is mounted in a three 
dimensional trapezoidal support mount with quarter-inch 
plexiglass walls on the top, bottom, and sides. The front 
and rear are open to allow easy insertion of the copper 
plate. The dimensions and shape of the nount are given in 
Pigur2 2.7 Ail dimensions are outside wall neasurements. 

The flat copper antenna 1s slipped into the 
front of the plexiglass mount and protrudes through the 


back. A coaxial cable is connected to the back with an rf 
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Figure 2.7. TEM Horn Support Mount Disensions 


end-connector and hand-tightened scraw. Ses Basie 2 emtor 


a diagram of the assembly of the connection point. 


antenna is physically supported within ‘he mount by grooves 
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cut into and running along the length of each side. The 
grooves serve to fix the flare angl2 of the antenna at 
approximately 24.19. 

The dimensions of the mount were chosen for 
reasons of practicality, being the most appropriate for the 
material available and providing a reasonable flare angle 
and antenna characteristic impedance for signal reception. 
The mount also provides stabilty to the Horn under vaxrious 
weather conditions. Tina's , by placing the horn in the 
plexiglass supvort mount, physical and electrical 
Seg@eaectesistics are f1xed, assisting in correlating run to 
run measurements. The TEM Horn as it would be mounted on 
the imaging plane is irawn in Figure 2.8. The support mount 
Bomso@tesnown for clarity. 

The addition of the plexiglass mount should have 
negligible effect on the received signai. Figure 2.9 shows 
the response of the supported horn antenna to a pulse 
generated by the source. This compares favorably with the 
pulses obt ained by Hammond without the plexiglass 
mount. [{Ref. 11] 

As noted in Figure 2.8, the electric field at 
the aperture of the horn, assuming a uniform TEM wave being 


excited between the horn and the ground plane, is given by: 





Figure 2.8. 


2 Vo -jske 
E-:Te °| 4h 
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(a) Front View 


(b) Side View 


TEM Hocn Receiving Antenna 
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is O 
E=—~ e - (2.1) 


Where: Vo = amplitude of the incident voltage at 
tne Lreed pont . 


h = height of aperture 
K = wave number = 2TT/» 
a = dastance of observer fron the 


referance point. 
It can then be easily shown that the receiving 


antenna impulse response in the far field is: [Ref. 12] 


o> (G2) re e, WZ, § (t - r,/c) (2.2) 
wheres: § = impulse 
Z, = load impedance 


and the impulse response of the transmitting antenna is: 





iC ME 0 
Ce =-- - at 6 (t-r /c) 
where: rc = distance from the aperture to an 


observer on the image plane in the 
DOEesight direction. 


Thus the transmitting impuls2 response is simply 
the time derivative of the receiving impulse response. The 
impulse response is a delayed impulse, delayed in tine by 

EP and multiplied by a negative constant which is 
dependant on the width of the antenna aperture, Ww, and the 


load impedance Z,. 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
{ 6-26-81 N/A HQ 777 GENERATOR WAVEFORM 
MAXIMUM PEAK VALUE CVOLTS) 9.2528 
MINIMUM PEAK VALUE CUOLTS> “9.0932 
RMS VALUE <UVOLTS) 8.8796 
MEAN VALUE CVOLTS) 9.8808 


NUNBER OF WAVEFORMS AVERAGED # | 


“=F hOB OOS eeQ_OeOQeereseaez ee @ BS Cees 


E-1V0LT 


OPTIMIZATION VALUE = 6 


DIRECT WAVEFORN 





E-9 SEC 


8.58 1.08 1.58 2.80 2.58 3.08 3.89 4.88 4.58 5.98 
Pigure 2.9. 


Generally, the 
applicable eS. the Ie 
measurements: [Ref. 13] 

1. The antenna is not 
transmitting antenna. 
2. The flare angla of the 
recelve energy fron 
transmitting antenna. 
3. Wavelengths at 


a true TEM mode, 


a@mstortion. 


high frequencies 


resulting 


Typical TEM Horn Received Waveforn 


following limitations are 
Horn during nost TDRL 
Pie tue fap field of the 


TEM Horn is not wide enough to 


all positions along the 


may not propagate in 


in attenuation and 


39 





The net cesult is a dé@viation from the ideal 
results of the impulse responses as stated above. When 
Bequre 2.9 iS compared with Figure 2.4, as cOLte. on Of sane 
pulse as received by the antenna is noted. However, for 
TDRL purposes, the reproduced waveform is still reasonably 
faithful to the original source pulse. 

Cc. Receiving Antenna -- Monopoles 

Monopole antennas are also used onthe fTDRL 
scattering range. One of eleven monopoles, ranging in 
length from 85 millimeters to 285 millimeters, in 20 
Millimeter lengths, can be chosen. Each monopole is 
constructed of brass rod and is threaded on one end to allow 
for ease of connection to the receiving elements. 

The monopoles may be considered center-fed 
dipoles with their images in the ground plane. In effect, 
they are oscillating electrical dipoles in free space. 
Oscillation at frequencies related ae) the 
quarter-wavelengths of the monopoles is to be expected. 

Fagmres 2.10, 2.11 and 2.12 confirm the above. 
The effects of radiation jiamping, dispersion anda reflection 
due to the impedance discontinuity of the antenna are 
clearly demonstrated. From Figure 2.4 the pulse width of 
the source pulse is seen to be about 1.2 nanoseconds. This 


corresponds to a wavelength of about 360 millimeters. 
40 








TARGET: RUN DATE OIST ANT TGT REMARKS 
4 6-29-81 1.27 2 777 ACQUIRE 85MM ANTENNA RESPONSE 


MAXIMUM PEAK YALUE CUOLTS) 8.0577 
MINIMUM PEAK VALUE CUOLTSD -2.9702 
RMS VALUE <UOLTSD Q. 8237 
MEAN YALUE (VOLTS) @.90080 
NUMBER OF WAVEFORMS AVERAGED = 1 OPTIMIZATION VALUE = @ 
E-2U0LT DIRECT WAVEFORM 


3.88 


8.08 


-5.98 





9.59 1.00 1.30 2.88 2.58 3.90 3.50 4.88 4.589 3.988 
Figure 2.10. 85mm Monopole Pulse Response 


Figure 2.10 is the short-pulse response of the 
85mm monopole. This antenna is almost exactly one 
quarter-wavelength in height in relation to the source 
pulsewidth. It correspondingly produces a nearly pure 
Sinusoidal variation, damped in tine, as predicted by 
Classical antenna theory. 

Figures 2.11 and 2.12 are the observed responses 
for the 185 millimeter and 285 millimeter monopoles. As can 
be seen, the antennas tend +o resonate at there 


quarter-wavelength periods of 2.5 and 3.8 nanoseconds 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
3 6-29-81 1.27M ? 777 ACQUIRE 185MM ANTENNA RESPONSE 


MAXIMUM PEAK VALUE CUOLTS) 9.8859 
MINIMUM PEAK VALUE CVOLTS) -8.8907 
RMS VALUE (VOLTS) 8.93298 
MEAN VALUE (VOLTS) 9.8800 
NUMBER OF WAVEFORMS AVERAGED = 1 OPTIMIZATION VALUE = 8 


DIRECT WAVEFORM 


E-2V0LT 
3.88 
8.938 
-5.98 
E-8 SEC 
8.90 9.20 9.40 8.68 8.8 1,00 
Pigure 2.11. 185mm Monopole Pulse Response 
respectively. Although the amplitudes of *hese two 


monopoles do not differ greatly (less than .08% between 
their maximum peak values), the longer antenna tends to 
distort the incident wave to a greater extent. This is most 
readily observed when compating the root-mean-square 
yoltages of the two waveforms. The 185 millimeter antenna 
oprovides more effective power to the receiving circuits than 
does either of the other two monopoles. The 295 @ilbageter 
antenna dissipates 4 considerable amount of energy in the 


long period oscillations; the 85 millimeter antenna 1s *0o 
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small to acquire the full power of the transmitting antenna. 
There nay then be a "best" monopole length for a particular 
source pulse. What this length might be and whether the 
monopoles are selective of targets according +o their 


dimensions will be discussed more fully in Chapter V. 


TARGET: RUN DATE OIST ANT TGT REMARKS 
2 6-29-81 1.27M 12 777 ACQUIRE 285MM ANTENNA RESPONSE 
MAXIMUM PEAK VALUE “VOLTS? 9.8861 
MINIMUM PEAK VALUE CVOLTS) ~8.8934 
RMS VALUE (VOLTS) 9.8289 
MEAN VALUE CUVOLTS? 8.0008 
NUMBER OF WAVEFORMS AVERAGED = } OPTIMIZATION VALUE = 8 


DIRECT WAVEFORM 
E-2V0LT 


53.88 


8.80 


-3. 88 


E=8) SEC 
8.98 8.58 1.88 1.58 2.88 


Figure 2.12. 285mm Monopole Pulse Response 


All receiving antennas are coupled to the Signal 
Processing Group through a coaxial end-coupling fitted 
through the imaging plane. The connection point is 


diagrammed in cut-away in Figure 2.13. The antenna is 
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Figure 2.13. Receiving Antenna Connection Point 
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coupled as close to the ground plane as possible. Tits as 
Particularly important for the TEM Horn to prevent 
undesirable oscillations. As drawn, the threaded portion of 
the nonopole (or hand tightened screw of the TEM Horn) is 
jcined to a connector pin by butting the surfaces of the 
antenna and pin together. A coaxial cable is then connected 
to the end connector and the rf signal is transferred to the 
meeea¥ing circuits. 

This accangement has worked well. No 
reflections or excessive attenuation of the incoming signal 
due to the joint were observed. Gserliations on the TEM 
Horn were minimal. 

Table I lists the antennas currently available 
to the TDRL user. The significant dimensions and pertinent 
characteristics are reviewed. 

d. Targets 

The targets used on the TDRL scattering range 
are of various dimensions. Cuerenc.y, eleven are used in 
imaging experiments. They are all axisymmetric, closed 
surfaces with edges and vertices. The targets may be broken 
into three catergories according to general shape: 

io ieali-cyivaders 


2. Half-spheres 
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TABLE I 


TORL Antennas 





























ANTENNA DesGrre lr LOW | BEYENSIONS | COMMENTS 
| NOMBER | p__3e 7 
a DE Horn | Figure 2.6 | Broadband 
| 2 : Dipole as 85mm Oscillating Dipole 
ao «| pipole. | 10cm 1. 
4 | Dipole 125mm 
co Ta5a0 : 
6 Dipole 165na { 
7 Dipole 7 18500 | 
8 Dipole Ries 205mm a 
9 | Dipole | 225mm 
10 Dipole ae i 
q2 Dipole =| Ssh a | ae 


J. faet-cones 

Table II lists the targets and provides general information 
about each. 

There are six half-cylinders used as targets. 
They are all machined from aluminum stock. The largest is 
12-inches long with a 3-inch radius. The smallest is 
3.06-inches long and 1.5-inches in radius. 

A single half-cone is also used asa target. 
The cone is machined from the same stock as the 


nalf-cylinders. 
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Four spheres, varying in radius from S-inches to 
1.5-inches are the four remaining targets. The spheres are 
basically styrofoam balls cut in half. Each half is covered 
with heavy-duty aluminum foil to provide the proper 
reflection of transient EM. 

Bach of the targets, when placed on the ground 
plane, appears to the incident pulse to be mirrored. 
Therefore, the receiving antenna sees not the half-shapes, 
but the full cylinder, cone cr sphere. 

No particular problems were noted when using the 
targets except for the need to anchor the spheres with 
weights during inclement weather. Otherwise, they have a 
tendancy to move on the plane due to the wind. 

e. Image Plane 

The final element of the Image Plane Group is 
the image plane itself. It is a square structure of welded 
Reynolds Type AN-190 aluminua sheet. Tt measures about 11 
meters on a side for an area of about 121 Sguare meters. It 
PemcOur mal limeters thick. 

Based on the above, the range clear time, which 
establishes the low frequency cut-off of the measurement, is 
jreater than 36 nanoseconds. The clear tine establishes 
the first arrival at the measurement point of a reflection 


from objects surrounding the plane and the end of the plane. 
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TABLE II 


TDRL Targets 


ee nen me he ee a fe ne ee ee 
| TARGET amo t Gee =a at MATERIAL | 
1 at 37-4 ~~ Aluminum Stock 
Cylinder 5" | 1p meter ae Aluminum Stock | 







2 
3 Cylinder 6" Aluminun StOCK 
= 4 { Cylinder Cx Ae! 6.06" | Aluminum Stock 
a eee rae aera | 
5 | Cylinder | as 3,00" | Aluminum Stock 
6 {| Cylinder es esis 3.06" Aluminum Stock 
{ 7 Sphere Sy Al and Styrofoan 
al 8 Sphere | yt Al and Styrofoan {| 
Es i —— a ed 
be 3, Sphere i 3" | | Ail and Styrofoan | 
SSS SSS pS SS SS SSS SS SS See 
4s 19 Sphere t ipo { { Al and Styrofoam {| 
eee ae oe ae ee ef | 
11 Cone ae vo. 5.94" Aluminum Stock { 
pe as ee eteereneenenseeamenioaeneinad 


The main purpose of the image plane is to 
isolate the targets and the antennas from the underlying 
structure of ene SDR. Instrumentation cables are 
introduced from under the ground plane and are thus kept 
from interfering with the measurement. The plane also 
Serves as the primary support structure for both the 
transmitting and receiving antenna. 

No significant problems were encountered with 
the image plane. Jowever, som2 disadvantages of this plane 


arrangement are that only objects having a symmetry plane 
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can be measuted, and polarization and incident angles are 
limited. But, for the work intended, the image plane serves 
its purpose well. 

3. Sampling Receiver and Signal Processing Group 

The Signal Processing Group consists of three 
elements: 

1. Digital Processing Oscilloscope (DP0O) 
2. Mini-Computer 
Saeueard Copier 

The Signal Processing circuits are connected from 
the output of the receiving antenna to the input of the DPO 
Via 2.58 meters of RG-8A/U coaxial cable. The cable 
provides some attenuation and delay of the source signal 
provided by the impulse generator. Additionally, a trigger 
Signal is sent to the DPD from the pulse generator via .54 
meters of RG-53C/U cable. 

The lengths of the signal and trigger cable are 
chosen for convenience and for minimum attenuation. They 
May be varied for special purpose applications such as to 
accommodate a different trigger source by providing a 
necessary delay. Care must pe taxen not to introduce too 
much delay or the received signal will not be windowed on 


the CRT display of the DPO. The maximum practical length of 
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Gadling, including that to the transmitting antenna from the 
pulse generator, should be less than 20 meters for the 
miherent equipment configuration of the TDRL. 
a. Digital Processing Oscilloscope 

The Tektronix Digital Processing Oscilloscope is 
a computer compatible, analog to digital systen. Less 
comprised of two elemental units; a Tektronix 7704A General 
Purpose Oscilloscopes Syst2n, and a P7001 Processor. 
Modification of the 7704A Oscilloscope by the addition of 
various plug-in modules allows for tailoring of the whole 
DPO system t9 the specific needs of the IDRL. A cCULSOLy 
description cf the various elements of the DPO, as 
applicable to the TDRL, follows. 

The 7J704A Oscilloscope system is composed of two 
Marts; the 07704 Display Unit and the A7704 acquisition 
wot. The Display unit provides results of signal 
acquisition, both real time and as stored in the P7001. ie 
does this by a visual interface with the DPO in the form of 
an Oscilloscope. cadtteomauly, the CRE readout 2s cazrice 
over into the DPO to provide appropriate scaling on the 
displays. 

Trem ter Ol acguacit1on Unit, Sthrowgh its. @odulac 


plug-ins, provides the means t9 acquire the analog signal. 
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The plug-in units used in the TDRL are the S-6 Sampling Head 
and the $-53 Trigger Recognizer, INSEEted into 42 Psi2 
fpR/Sampler. 

The 7512 provides for measurement of recurring 
fast-rise tine signals; i.32., the source pulse and the 
target back-scatter. It features the means of determining 
vertical deflection factors in millivolts; horizontal sweep 
factors in seconds; a HIGH RESOLUTION switch which reduces 
the waveform nois® and jitter by signal averaging (10 
averages per sampl2 point); a TIME-DISTANCE scal@ which 
allows one-way distance measurements in air dielectrics; and 
a locate switch which increases the tine/division and 
intensifies a portion of the display to locate the tine 
window relative to the total waveforn. Pemaliy, che 7S 12 
provides the means for the inserting and interfacing of the 
plug-in heads into the DPO. 

The S-6 Sampling Head is a 50-ohm loop through 
input sampling unit. It provides the means of inputting the 
Seaiaeeo. tacerest ¢95 the DPO. This unit has a specified 


rise time of 30 picoseconds or less and a bandwidth 


equivalent to de to 11.5 GHz at 3 dB down. Tt has a Signal 
Bemge Of +1volt to -Ivolt (dc to ac peak) or tvolt 
peak-to-peak. The S-6 establishes the baseline measuring 


capabilities of the IDRL. 
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me =S-53 Trigger Recognizer vermits use or the 
7812 aS a general purpose sampler. It produces a stable 
trigger signal from input signals from DC to 1 GHz. eS 
Capable of recognizing signals having amplitudes ranging 
from 10 millivolts ts 2 volts peak-to-peak into 50 onms. A 
trigger pulse from the pulse generator is fed into this unit 
at the same rate as the source signal. Irigger to signal 
delay is specified as 15 nanoseconds or less. 

The second element of the DPpoO, tne P7001 
PLOCeSSOL, provides the lee FLO the DPO to the 
mini-computer. This unit receives th2 analog signal 
acquired from the Acquisition unit and digitizes and stores 
or discards the information, along wit appropriate scale 
Factors, as selected by the operator. The method used by 
the digitizer is pseudo-randon sampling with one sample 
taken every 6.5 microseconds. This is the maximum rate of 
sampling. Any transient longer than 5 milliseconds or 
COHERENTLY REPETITIVE signal up to the frequency limit 
specified by the Sampling Head that can be displayed on the 
D7704 Display Unit can be stored in internal memory along 
with its scale factors. The information stored can then be 
redisplayed on the CRT and/or sent to the minicomputer for 


processing. There are a maximum of 512 samples taken 
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horzontaily (99-bit word) over waveforn. Although either X- 
or T- rsp May be specitied for this axis, the TDRL 
primarily employs the time axis. Vertical resolution of the 
Processor is 1024 levels (10-bit word). This axis is used 
to measure the voltage of the back-scattered signal. 
b. Tektronix 4052 Mini-Computer 

The minicomputer used in the TDRL is a Tektronix 
4052. It is a high performance, integrated graphics systen 
employing LSI, bi-polar, 16-bit technology. The unit has a 
64k-byte memory with a 300k-byte magnetic cartridge built 
iA). Extended BASIC is the high level programming language 
employed to perforn orocessing functions. Although this 
unit has many features worthy of note, the above are the 
Salient ones for the TORL. 

SeeHerd Copier 

The Hard Copier used is the Tektronix 4631 Hard 
Copy Unit. If provides permanent, dry copies of the graphic 
and alphanumeric information displayed on the CRT storage 
screen of the 4052 mini-computer. 

All units described for the Signal Processing 
Group were especially selected to perform as an integrated 
system requiring minimum interfacing. The GPIB bus is the 
Boarmiseon transferring information from one unit to the other 


in the TDRL. 
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This chapter is devoted to the general derivation of the 
equations which are applicable to the «time-domain analysis 
of the impulse and ramp responses of axisymmetric, metallic 
bodies which are illuminated by an axially directed incident 
electromagnetic field. “here sources providing detailed 
derivations are numerous and/or adequate, the derivations in 
this chapter will be general, providing only the major 
results of interest with specific constraints noted. 

The time-domain integral equations which provide 
solutions for the transient electromagnetic problem wiil he 
considered first. They will be followed by a discussion of 
the measurement equations and procedures employed in the 
Dont. Finally, various methods used to reduce the noise 


error in acquired signals will be discussed. 


A. TIME-DOMAIN INTEGRAL EQUATIONS 

One primary problem in transient analysis of target 
impulse responses is to determine the shape of an unknown 
Object when the incident fieid and the response of the 
Bee@eerer to that field are a priori information. Numerous 


time-domain technizues have been developed which perform 
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mms LSunction. As noted in Chapter I, tha earliest method 
used involved the physical optics response of a body to an 
Meecrzdent electromagnetic field. The surface currents 
produced an approximate impulse response that was 
essentially the second derivative of the projected area 
function of the scatterer. [Ref. 5j] However, this approach 
was limited in its accuracy in that the relationship between 
the impulse response and the derivatives is exact only at 
the leading edge of the scattered field response, a single 
point in time. Interactions between currents on the target 
pody which continue to radiate for a significant time after 
the incident field leading edge has propagated further in 
Space alters the backscattered tesponse of the target. 
Therefore, for an accurate solution to the problem, these 
additional currents must be of consideration. 

Bennett {Ref. 14] first proposed a direct method which 
"corrected" the response by computing the currents flowing 
on the scatterer surface. From these currents, the 
Bea@eered field could be calsulated "exactly". 

In general, the remainder of this section follows the 
format of the method used by Bennett in the derivation of 


integral equations for exact solutions. 


ahs 








Pee expression for the wagnetic field 4 at an 
arbitrary point in space (not on the scatterer surface) is 


given by: 


x = eae ~ = ’ -~ 
H(r,t) = qr) + 7 f (I< + a) J36G a) ap ds (oe) 


a 


a 
wheres H(r,t) 


Gotal smagnet ie@imeld@ae (r,t) 





a' (Z, t) = incident magnetic field at (r,t) 
J(c', +t) = surface current at (f',¢) 
= o e 
ig = position vector EO) the 
observation point 
ee = position vector to the 
integration point 
R = |E-E'| 
cab 
A — rer! 
“Tae 
—~ ¢ = time ig light-meters (one 
light-meter is the time it takes 
light to travel one meter) 
Y = t-R = retarded tine 


By specializing the arbitrary space point to T, a 
point on the surface of the scatterer, and then applying the 


Boundary conditions (1i.e., by causuality, incident field is 
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@etoePETOr to arrival of the incident pulse ata sample 
point), an integral equation for the current density on the 


Surface of the scatterer can be round: | 


ee ws ae 5 eee 1 2 1. Apes. toa 
Crs €) ~ 2a. x oH (Se) iz a f a x i a R al Cree 
cuanl eds! Ce) 
where: Tt = position vector to a point located on 


the surface of the scatterer 


URese Vector NOrEMal tortne suctrace 


3? 
it 


The first termon the right-hand side of equation 
3. 2) is the source term. It represents the direct 
mmehucnce’ 9f the incident field on the current at the 
observation point (E,t). iG is the physical optics 
approximation for the surface current when applied to the 
iliuminated side of the scatterer. The integral on the 
ergne side of equation (3.2) represents the influence of 
Si@erents at other surface points on the current at (£,t). 
Note that the influence of other currents on the current at 
(r,t) 1s delayed by R. This allows the surface current 
Gensity equation +9 be solved by an iterative numerical 
procedure in the time-domain, rather than the familiar 
frequency-domain matrix inversion process. 

Once the surface currents have been found, *he 


far-scattered field can be calculated by using: 


a7 





S,> 1 -) ze = —— re 
rH ye = 7— ae ey ee) X at as! (So 


where: aN dastance to the far—-fiela scbhserver 


er = unit vector from the integration point 
to the far-field observer 

W2oh the Substitution of the Surface current 

expression into the above equation, andthe assumption that 

the incident field is a ramp, the result for the backscatter 


Jeerection is: 


r HS rot ap = re: : = ~ 
Pe p 6 ) TES S(t. Ja, ae = Vin (Jon (r,t) x a} ds! er) 


ly 
where: fC i. backscatter ramp response of the 


oO 
target 
S(t.) = physical optics silhouette area of 
Ene wanget 
Bo = distance Of the far-field observer 
from the origin 
= + 
is ag aS 
- 2 
: |H™ | 
aN =J,, correction currents resulting 


from the incident tramp waveform. 
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This result provides the exact relationship between 
the target response and target geometry. A particular 
unknown Sample value at a given point in space and tine is 
determined by the exciting field at that same space-time 
point and by the scattered fields from earlier, more distant 
focatlons. Note that the interactions between unxknown 
Samples are displaced in time by an amount equal to the time 
reguired for a field to propagate between the samples at the 
speed of light. The unknown Samples can be solved at any 


time step, provided all sample values at earlier times are 


already known. Thus, by determining the unknown samples, 
Ds Cs, the GQOrrect21on ‘currents ah: and adding the 


Semtcibution of the target area function, S(t.) the target 
ramp response can be utilized to recover information 
involving the unknown dimensions. 
2. Inverse Scattering Solution 

The targets employed siya the TDRL are all 
cotationally symmetric scatterers, Similar to the one 
diagramed in Pigure 3.1. The general constraints applicable 
to the problem are that the scatterer is symmetric about the 


Z-axis, the incident field is axially incident, and the 


far-field is computed in the backscatter direction. 


29 








Pigure 3.1. Geometry of Rotationally Symmetric Scattering 
Probles 
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feremermeOur Df a2 rotationally symmetric target is 
completely described by the radius vector e (2) which 
varies as a function of 2z. Bic moLOVvecceararce  fUNGCT1ON Ln 
the plane orthogonal to the direction of the propagation of 
the field (z-axis) can be expressed simply as: 

5 = 197 (2) 

Substituting the projected area function into the 
aquation for tne backscatter ramp response of the target and 
solving f° 5 (Zz), an inversion equation for the 


rotationallly symmetric case is obtained: 


1 WZ 


ac seg 3 2 ' t ae 
o(2) = (2rgg(F.t) - grag § egltht) x aj} ds] say) 


i (S25) 


Thus, O (z) is given in terms of the measured ramp response 
(a field measurement), and in terms of correction currents 
at earlier times, which have been previously computed or are 
Known to be zero. 

In the TDRL, an iterative approach is employed using 
astimates of the entire target geometry to solve for the 
eonmcoure fUmc<1on. [The procedure involves five steps: 

1. Neglecting the surface integral providing the 
correction currents, obtain the physical optics 


estimate of a 





2. Determine the correction currents using the numerical 
solution of equation (3.2) £0r sSuEzteace scurrent 
density. 

So APpieyecne §«SCOrrection currents to equation (3.5) to 


obtain a next estimate of QP (2). 


e,(2) = [2r) He (F,2))1/? 


4. Compare the new value e (Z) @2eh P (2) EOS ees 
2 

the change is less than some small number, an error 

Factors The normalized mean sguare error is defined 


ass 


(H3(t) - HS (4))? 
k-1 


R = 
€ = a — sree 
K S Z 
(Hp (t) ) 
where: a = the value to be minimized 


in order to obtain the 
nost accurate shape 
description 


the kth sstimation of the 


re 
, a) 
ct 

tt 


backscattered field. 
Thus, the difference between the astimated and actual 


backscattered fields are directly linked to the 
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difference between successive contour functions, and, 
therefore, successive shapes of the body. 
5. Repeat the iteration orocess until a minimum error is 
obtained. 
The above procedure was used successfully by Morag 
(Ref. 15] in the development of an algorithm involving 
theoretical analysis of targets. This algorithm has been 
integrated with the overall IDRL algorithm for inverse 


scattering. 


B. BACKSCAITERED WAVEFPORA MEASUREMENTS 

The TDRL algorithm is designed to minimize the effects 
of the impulse responses of the transmitting and receiving 
antennas. Figure 3.2 diagrams the parameters involved. 

As can de seen, there are two measurement paths of 
interest. The first is the Direct Path between the 
transmitting and receiving antenna. Micon ys taceors 
affecting the transmitted vywaveform are the intervening 
distance andthe impulse responses of the receiving and 
transmitting antennas to the source pulse. The second path 
is the Augmented Path, from the transmitting antenna *o the 
target and back to receiving antenna. Factors effecting the 


transmitted pulse in this path include the antenna impulse 
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IMPULSE WIRE 


FREE 
GENERATOR ANTENNA SPACE DIRECT 
_ |DISTANCE D PATH 


(S(w)) cHo (w)) 





CH, ( w)) 








FREE _ FREE 
SPACE TARGET SPACE 
DISTANCE D Z DISTANCE 2D 
(H (w)) 
CH, (w)) CH, (Cw) ) 
: RECEIVING 
ANTENNA 
(H™ (w)) 
NOTE: 


Effects of wire transmitting antenna 


reflections at receiving antenna considered 
small -- neglected. 


Pigure 3.2. Signal Parameters for TDRL Measurements 


64 





responses, the impulse response of the target and the 
corresponding intervening distances. 

Three waveforms are measured: the Direct Waveform, the 
Incident Waveform and the Augmented Waveforn. During the 
measurements, h' (t) and h* (4), responses due to the delays 
of the intervening distances, are windowed out by using the 
time-range gating feature of the DPO. 

The Direct Waveform 1s measured first. The transmitting 
and receiving antennas are boresighted. No target is on the 
imaging plane. The measured wave is then the convolved 
responses of the source, and transmitting and receiving 
antenna impulse responses. By using the properties of 
convolution, tne frequency domain representation converts 
the convolution process into a product. The frequency 


domain representation cf the Direct Waveform is: 


Aced®) = s(ed*®)Ht ed”) HT (21%) 


wheres Rien ae Direct Waveform 


s(e°’) = source pulse 

Hye") = transmitting antenna response 
we’, SW _ tas 

H(e"" ) = receiving antenna response 


The effects of the intervening distance have been time 
gated out. Figure 3.3 is a representative Direct Waveforn 


measured in the TDRL. 





TARGET: a ee aan OIST ANT i REMARKS 


Mest l.27t 1 CYLINDER ACQUISITION AVERAGE 
MAAPMUM PEsmme URIBUE. sce cc cence +289.21 mv 
MINIMUM PEAK VALUE. cccccccccce “96.55 av 
RPS OMIMUEs cs ccccsoeccccccscrece +98.58 mV 
MEAN VALUE. . cece secre rccrvecee 8.80 mV 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = @.35 





i i a 
YG Ss 


8.58 1.88 1.58 2.08 2.50 3.88 3.59 4.08 4.38 5.98 nS 


Figure 3.3. Typical Direct Wavefora--Cylinder 
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The Incident Waveform 1S measured next. The imaging 
Blane das as set-up f9r the Direct -‘Waveforn, but the 
time-gating is further delayed by dat to eliminate the 
effects of the distance intervening between the receiving 
antenna and the target site. The resulting measured 
waveform is the Direct Waveform delayed by dt: 


B(eJ) = Ate”) o (jwdAt) 


wheres; Bie” ) = Incident Waveform 

rigure 3.4 provides a representative time-domain Incident 
daveftorm measured in the TDRL. 

The final waveform measured 1s the Augmented Waveforn. 
It is a convolution of the Direct Waveform with the impulse 
response of the target and Summed with the Incident 
Wavertorn. For this measurement, the transmitting and 
receiving antenna remain boresighted. A target is vlaced on 
the image nvlane. The measured frequency domain 


representation is: 


Cel”) = [B(eJ%) + aceJ%) HS(ed) 1 _Juddt (sae) 
wheres Giana) = Augmented Waveforn 
H (e@°) = Target Response 


Figure 3.5 provides 2 representative time-domain Augmented 
Waveforan measured in the PDs L for a axisymmetric 


half-cylinder. 





TARGET: oe DATE DIST ANT TGT REMARKS 
Geral lcerm 1 t CYLINDER ACQUISITION AVERAGE 


MAXIMUM PEAK URLUE.. pc eee nceee 4.40 mU 
MININUM PEAK VALUE... cc ecenee -3.83 mU 
RMS VALUE. ccceccccnccrcrvcsses #1.92 mV 
NEAN VALUE... cc ccscccscccceaes 0.00 mV 
NUMBER OF WAVEFORNS AVERAGED = 21 OPTIMIZATION VALUE = @.5 


INCIDENT WAVEFORM 


Vv SSR REE 
“SEP 400 
; ~ 5340 Ree 
PIN AT TET A 
~ REECE 


@.58 1.80 1.50 2.08 2.50 3.989 3.58 4.398 4.58 5.880 nS 





4,20 








2.80 


8.88 





-2.088 


Figure 3.4. Typical Incident Wavefora--Cylinder 





TRRGET: RUN DATE DIST ANT TGT REMARKS 
isin Cle aieari tf i CYLINDER ACQUISITION AVERAGE 


MAXIMUM PEAK MNGREUIE Gece e606 6.0 6 0 6 4.93 mV 


MINIMUM PEAK VALUE... cccccee “6.76 mV 
Tt OMNES aisle + sca closes es eee cces 42.76 mV 
MEAN VALUE. . cece reer c cre cveane 9.08 nV 
NUMBER OF WAVEFORMS AVERAGED = 23 OPTIMIZATION VALUE = 9.5 


9.80 





Q.58 1.388 1.50 2.808 2.58 3.89 3.58 4.98 4.589 5.88 ns 


Pigure 3.5. Typical Augmented Waveformu--Cylinder 
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AS can be seen from eguation (3.6), the target response 
is not independant of the effects of the impulse responses 
of ‘the antennas. To eliminate this dependance, the 
backscattered waveforn is determined and then the remaining 
antenna terms are theoretically eliminatei by a division 
process. 

The Backscattered Waveform is the difference of the 


Augmented and the Incident Waveforms: 
p(o2%) = ce) - B(eI%) 
_ jw, WS -,Jw 
= A(e- ) H™ (e-) 
eeaure, 3.6 is the time-domain Backscattered waveform 
resulting from the subtraction of the waveforms in Figures 
Ba Jeand 3.4 respectively. 


By next dividing the Backscattered Wave by the Direct 


Wave, all antenna influences are eliminated: 


jw 
H?(eJ%) = ame 
mice.) 
Toe target impulse resvonse, h> (+), may now be determined by 


deconvolution of the target impulse frequency domain values. 
Figure 3.7 is a representative time-domain target impulse 
response for a cylinder. Noise errors have been minimized 


Dy an adaptive filtering method +o be described shortly. 


70 





TARGET: RUN OATE DIST ANT TGT REMARKS 
Preeti -ss. laerm 3% 1 CYLINDER ACQUISITION AVERAGE 


MAXIMUM PERK VALUE... ccc cceceee 44,65 mV 
MINIMUM PEAK VALUE... ccveces -4.44 mV 
Rita VUE ss ccc secs ccs cccescces 1.71 mV 
MERN VALUE... ccc cece scree cccccs 8.98 mV 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8.5 


BACK SCATTER WAVE FORM 


a 
CCAR ECCT 
OC a 
PCI A 
CONCEP 
a 


6.58 1.08 1.580 2.909 2.58 3,30 3.58 4.98 4.58 3.88 ns 








Pigure 3.6. Typical Backscattered Waveform--Cylinder 
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TARGET: RUN DATE OIST ANT TGT RENARKS 
esse oleeisearh 1 i TYPICAL CYLINDER RESPONSES 


MAX INUM PEAK UALUE. @®eeeeseeeseesteoe Q. 8804 hy, 
MINIMUM PEAK VALUE... .cecceer. “8.0004 mV 
RAS VURLUE. @@eeeeesede#es+rieseee#nw1eeteee#ee#s 0. veG2 mV 


MEAN VALUE. eee 8.08000 MU 
NUNBER OF WAVEFORMS AVERAGED = l UPTIMIZATION VALUE 2 18 


SEO SOO O@OB* OO SBT SBeG* ee OQ GaSG ae 22 SS SB SSB SO FG SF Sf SF SH S SH SH Sf SF S SH SP SF SB SL SF VSHVPFPQS#@ QDS SG = 


8.9804 


8.8882 


8.8888 


-O. 8682 





@.58 1.08 1.58 2.88 2.58 3.88 3.58 4.80 4.50 5.88 nS 


Fugure 3.7. Ly ical Ranges: Impulse Response With Noise 
ritering=-Cylinder 


72 





C. NOISE CONSIDERATONS 

As noted in the previous section, deconvolution of the 
target backscattered impulse response, less the iapulse 
responses of the antennas, is a central operation to the 
effective use of the IDRL. The reliability of the acquired 
signal to be a true response of a target to an incident 
field bears directly on the accuracy of the final results as 
determined by solution of the inverse integral equations for 
time-domain analysis. Any errors in the data representing 
the Incident Waveform and Augmented Waveform degrade the 
cancellation of the division process and introduce errors in 
the backscattered target response. Also present at the 
input to the DPO is noise due to thermal electron excitation 
in the transmission Lines and antenna surfaces as well a 
that due to external atmospheric, cosmic and man-made noise 
sources being received by the horn antenna. When the 
backscattered impulse response is transformed to the 
time-domain, high frequency noise in most cases forms the 
dominant feature of the waveform. An example of a typical 
impulse response of a cylinder is given in Figure 3.8. 
Noise errors have not been removed. Figure 3.9 is the raap 


response for the same target. 


73 





TARGET: RUN DATE OIST ANT TGT REMARKS 
feeeeie ss ilearm i 1 TYPICAL CYLINDER RESPONSES 


Meanwie: PEAR VALUE. ...c.ccceee “O-0I24 mV 
MININHUM PERK UP US ra “0.0314 my 
RMS WSL USS AA 8.0126 mV 


MEAN VALUE. eooe 98.9008 mV 
NUMBER OF WAVEFORNS AVERAGED = | OPTIMIZATION VALUE = 2 


<> =D 2p GP 2D aD aD SD Oe oD 4D OD o> SD OD OD OD 6 Oe SD SD ED ED SE SD ED ED ED DS ow OD ES SO Se S & & SED o® 4D o> SD SED SD OE SD SD OD OD o> OD om OD 4 OO ee ee ee ee Se eS ee ee ee oe 























9.9209 HH HL oA YATE ae | \ Hn 
2.2009 julie ll! til i Mt Hi " vil 


a | 


9.59 1.90 1.50 2.88 2.58 3.98 3.50 4.80 4.50 5.890 nS 





Pigure 3.8. Impulse Response--Noisy Signal 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
Peece=iuee, leer 1 1 TYPICAL CYLINDER RESPONSES 


MASIMUM PERK VALUE... cccscncee 4353.66 mU-M 
MINIMUM PEAK VALUE... cece eee ee “700.90 AU-M 
RMS VALUE. cece ccccccccccvecccs§ P56 se IS MVM 
MERM URL US ccc ccc ees ccrccecess§ @2Sd.42 MUoM 
NUMBER OF WAVEFORHS AVERAGED = | OPTIMIZATION VALUE = 8 


De ee we SBS SSP 28 SOE SSS SES S|] SBS 2] SO 8 SF SP BPP ESTES SSS SCS OCS | SSF SP @ SF 2S SE 2S BPS Be SS 4S SS 2S - | =P = =P. = =! ab os = 


RAMP RESPONSE -- TIME 


LLL 
PCC 
ACL 
OA 
PCC Le 


0.58 1.80 1.50 2.80 2.58 3.38 3.50 4.08 4.58 5.88 nS 







Figure 3.9. Ramp Response--Noisy Signal 


ies. 





Mmiecee leet, ctyYpical sources of errors in the Incident 
and Augmented Waveforms are sampling errors and noise in 
Signal acguisition, leakage and aliasing errors in the 
Fourier Transfomations, and rounding errors in the computer. 
The error of primary concern is that involving noise in 
signal acqums Lt 270n and in particular +he enhanced 
high-frequency noise resulting from deconvolution of the 
antenna responses t9 obtain the impulse and tramp responses 
of the target. This noise has been directly traced to the 
Structure of the data within the input array, as will be 
noted in Chapter IV. 

Two methods are employed to reduce noise to acceptable 
levels: acquisition averaging and an optimal compensation 


technigue. { Ref. 16] 


Acquisition averaging is aimed at reducing signal 
acquisition errors and noise by averaging several 
acquisitions taken within a short period of time. During 
the process, a more serious error may be generated, that due 
Powe. ghal Gmiit with time. Signal averaging has been most 
useful in ensuring a good target ramp response. Tt has had 
noticeable affect on making the information content of the 


impulse response more available. Pigure 3.10 demonstrates 


76 








an impulse waveform with 210 acquisitions for the same 
maede=. aS for Figure 3.8. Theoretically, the RMS noise 
level should be reduced by the square root of the number of 
acquisitions, N. N = 10 in Figure 3.8, and N= 210 in 
Figure 3.10. The noise reduction should be about 4.5 times. 
Note that the noise level in Pigure 3.10 is indeed about one 
Fourth that in 3.8. 

Figure 3.11 is the corresponding ramp response. 
Signal acquisition only is employed. 

2. Optimal Compensation Technique 

The optimal compensation technique involves’ the 
design of a compensator (deconvolution) function operating 
on the convolution output of the Incident and Backscattered 
Waverorms, to produce a relatively noise-free estimated 
Impulse Response for the target. The compensator is applied 
to the convolution product of the target impulse response 
and the Direct Waveform to produce the required 
deconvolution result. The design as used in the fTDRL 
involves an iteration process on a Single variable with a 
Snman-in-the-loop" to determine optimal values of the 
variable. 

Pigure 3.12 diagrams the design of a frequency 


domain optimal compensator. The compensation principal in 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
WP eseiieet lezen 1 i CYLINDER ACQUISITION AVERAGE 


MAXIMUM PEAK VALUE... cence eee 8.0094 mV 
MINIMUM PEAK VALUE... cee ceecee “O.0091 mY 
RMS VALUE. ccvccvcccccccvceseee 98,0033 AY 
MEAN VALUE... cee cer ecccveceves Be-9000 MV 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = @ 


8.0000 


-8.8859 





8.58 1.80 1.50 2.00 2.56 3.88 3.58 4.99 4.58 5.08 nS 


Pigure 3.10. Impulse Response: No Optimization--Acquistion 
Averaged 
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TARGET: RUN DATE DEST ANT TGT REMARKS 
we Seles leer 1 CYLINGER-ACQUISITION AUG 


MAXIMUM PEAK VALUE... --ccccece *82.97 mU-M 
MINIMUM PEAK VALUE. . ccc eeceees ~441.33 MUMM 
RMS VALUE... ccc cccrecrcccvecee 293.42 mMU-M 
MEAN VALUE... ccc ccc cere cvceee § TIPS. 46 MUNN 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 


Se Se eee ee Ge Ge @SS ee Ge eee SOS SSB SSeS SO PC SP SF SO GB eT HB eS KF SESS SF SSS SE 2S es eGses Sa @ SS @ @* &@ aes 





-209 





8.58 1.08 1.58 2.08 2.58 3.980 3.58 4.00 4.58 5.398 nS 


-469 





Pigure 3.11. Ramp Response: No Optigization--Acquistion Av- 
eraged 


ie) 





= 





=i 
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Figure 3.12. Frequency Domain Optimal Compensator Design 


y e ® * e ae 
deconvolution is to design a transfer function Q(e ) to be 


applied to the backscatter response Bie”) to yield an 


3S ow 


estimate of the target impulse response H, (e jeeeinat 1S 


Ho 
e 


(eJ%) = dled} Q(e7%) (3.7a) 


eee aces”) H? (22%) Q(e7%) (ety 


The main design criteria is to ninimize the energy 
in tne noise terms that introduce error in the estimate. The 


error energy, E is given by: 


2 S jw S 
OMe He (eo) - H- (e2") 1? au 


Te) (3:23) 
where: (Q.= frequency band of interest. 
: ae ee S$, ww : 
When the error energy 1S minimized, Heafe  j > will 


equal H> (2?) and the impulse responses will be identical, 
except for the addition of highly enhanced random noise. fn 


this case, the compensator transfer function is: 
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‘ee ao. | 
Q(e?”) = ia *. (eJ%) = HP (2) 4) (3.9) 


But, this will result in an infinite, unbounded value for 


0 


dw ; ; 
Q(e” ) as A(e ) goes to zero, ylelding a noisy form for 


fee?” ) . This noise can be limited by constraining Q(2 ) 
+o be a bounded function. Noting that the impulse response 
of the target is bounded by physical constraints, Lees 


possible to keep the energy due to the convolution of the 


impulse response and the designed transfer function bounded, 


Bees S 
2 jw jw 2 
BE. = 4 | He?*) Q(e2®) [% du S55) 
The problem now is to minimize & while keeping E 
eanite. The total energy in the backscattered wave can be 


defined as: 


E=E,+E., A> 0 (3.11a) 


ry} : : : : 
ft] HS (e3%) - HS (ed) [2 + a ace?) Q¢e7) | “Jaw (3.116) 


where: \= Optimization parameter. 
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Thus, by minimizing the total error energy through 
the selection of the optimal value of A, the error energy 
in the estimate is also minimized. The result yields an 
Optimum compensator oer (2) ror the unknown target impulse 
response ao aE The result is that the whole system output 
is low-pass filtered with mininal affect on waveform shape. 

Adhering to the design constraints listed above, it 


can be shown that the form of the compensating transfer 


function 9 (e”) for a given Direct waveform A(e>”) is: 
(ref. 17] 
ey eat (et (Ropes 


eee) 
[JACe?) | %+a] 


where; Ax (e>”) = complex conjugate 
From equation 3.7, the compensation deconvolution 


process yields: 


“o.4 pcet®)a* (eo?) SNS) 
=] 


[ A(eI”) |*+a] 


Figure 3.13 shows the effect that optimization has 
on the impulse response with no sgnal acquisition averaging. 


Omermm—aization is 0.5. The target impulse response, while 
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TARGET: a DATE DIST ANT TGT REMARKS 
salt=3i 0 i.276 1 ! TYPICAL CYLINDER RESPONSES 


MAXIMUN PEAK URLUE.....ceeeeee O-GO11 mY 
MINIMUN PEAK VALUE. @ee@eeeee7#*ese8e@ -0, 8012 mU 
RMS VALUE. eo@eeeseeeeseeeeoeeeeeepeeeoeeee8e 0. 8984 mY 


MEAN VALUE. eeeeveee 9.8000 mV 
NUMBER OF WAVEFORMS. AVERAGED = 1 OPTIMIZHTION VALUE = @.5 


IMPULSE RESPONSE 


3.0018 






8.8885 






Mit tat ah PE |p 
rete LULA Wa 
LV ae 
SU 


8.59 1.00 1.58 2.808 2.58 3.98 3.50 4.88 4.38 5.80 nS 


8.9888 


-3. 9805 





Pigure 3.13. Impulse Response: .5 Optimization--- 
No Acquisition Averaging 


ES 








still noisy, is clearly discernable. Figure 3.13 should be 
compared with Figure 3.11. 21gurs 3. shows tneser fect of 
optimization increased to 10. The noise is almost totally 
eliminated. The trade~off is that the optimization has 
reduced the amplitude of the response and spread it in 
frequency. This is a general observation. However, it 
appears that the degree of effect onthe waveform is a 
Pemerion™ of the ideal optimization factor. Chapter V will 
explore this concept more fully. However, it should be 
noted that if the optimization factor is large enough to 
reduce the large scale features of the target then it is too 
large and should be reduced--this is the "man in the loop's" 
OD . 

The combination of waveforn averaging and 
optimization on the impulse response is demonstrated in 
meaures 3.15 and 3.16. With a relatively small amount of 
optiwaization, the target impulse response is readily 
apparent. See Figure 3.16. Signal acquisition averaging 
also yields an important advantage. Since noise is reduced 
directly by the square root of the number of acquisitions, 
and because no filtering of any sort is employed for 
averaging alone, the detailed features of the actual signal 


are not eliminated as with filtering. The smaller amount of 
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Optimization needed results in Minimum effect by the 
compensator on the shape of the output waveforn. The 
trade-off is the increased liklihood of error introduced due 
to Signal time drift during the acquistion averaging period. 

Figures 3.17 to 3.21 show the effects of signal 
averaging and optimization on the Ramp Responses. Each 
figure is presented in the same order as were the Impulse 
Responses. Signal averaging is very beneficial in obtaining 
an accurate ramp response. Optimization has a greater 
tendency to effect the wave shape of the ramp response than 
the impulse response of the target. But Less optimization 
is reguired to reduce the error noise. It is generally 
possible to use the ramp response without optimization. The 
Main effect of optimization is to spread the frequency 


content of the ramp response. Amplitude effects are minor. 
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Figure 3.14%. Impulse Response: 10 Optimization--No Acguisi- 
tion Averaging 
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Pigure 3.15. Impulse Response: No Optimization--210 Acqui- 
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Pigure 3.16. Impulse Response: .5 Optimization--210 Acquis- 
tions Averaged 
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Figure 3.17. Ramp Response: No Optimization--No Acquisition 
Averaging 
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Figure 3.18. Ramp Response: .5 Optimization--No Acquisition 
Averaging 
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Figure 3.19. Ramp Response: 10 Optimization--No Acquisition 
Averaging 
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Pigure 3.21. Ramp Response: .5 Optimization--210 Acquisi- 
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IVY. OPERATING SYSTEM COMPUTER IMPLEMENTATIO 


A. INTRODUCTION 

A set of software programs called the "Operating System" 
controls the input, processing and output of information 
obtained in the TDRL. It is the quality and throughness of 
the Operating System that directly determines the quality of 
the results obtained. In tais chapter, there will bea 
discussion in detail of the iliaplementation of the Operating 


System. It will begin by discussing the physical structure 


of the the host processor. A brief overview of the 
Operating System Wo, LOL uow, defining the basic 
Sub-programs that make-up the total algorithn. Finally, a 


detailed description of each program will be provided, 
emphasizing the method of implementation of the numerical 
Solutions, where applicable, and giving details of a 


procedure to run the individual programs. 


B. THE TEKTRONIX 4052 GRAPHIC COMPUTING SYSTEM 

The host processor used in the TDRL is the Tektronics 
4052 Graphic Computing System. It is a versatile, stand 
alone microprocessor. AS such, its capabilities and 
limitations are important driving functions that shape the 


structure of the Overating System. 
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The 4052 has three memories. They are: 

1. Read Cnly Wemory (ROM)--a permanen* memory that 

contains the System's intelligences. 

2. Random Access Memory (RAM)--working memory of the 
processor's CP9. 

3. Line Buffer--small temporary memory that allows data 
to be written and edited on a display screen prior to 
releasing it to the RAM. An important editing 
feature. 

Built-in peripherals for the Graphic System (GS) area 
Smart keyboard (primary input), a CRT visual display screen 
(primary output) and a magnetic tape unit (mass storage). 
Three extenal peripheral devices extend the versatility of 
the system. They are a Tektronix 4631 Hard Copy Unit (makes 
paper copies of display information), a 8924 Digital 
Cartridge Tape Drive (additional mass storage) and a Digital 
Processing Oscilloscope (data acguisition). A block diagram 
for the system is shown in Figure 4.1. [Ref. 18] 

The Processor is the main computing device. ie 
corresponds to the Cantral Processing Unit (CPO) found in 
larger systems. It maintains the "firmware" that allows the 
GS to direct system operations, decode instructions and 


perform arithmetic and logic operations. 
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Figure 4.1. Block Diagram of TDRL Graphic System 


The General Purpose Interface provides the means for the 
system to communicate information to or receive data from 
any external peripheral that is compatible with IEEE 
Standard 488-1975. MeEsm nt ertace 1S a  ~bic-daralled, 
byte-serial system capable of transfering AScrIt cod2 or 
machine dependant binary code up to 8-bits wide at a maximum 


rate of 250-k bytes/second. 
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A magnetic tape data cartridge is used in the internal 
magnetic and peripheral tape drive units to provide up to 
300-K bEytes of additional mass storage per tape. Programs 
and data may be stored on tape in specifically defined 
"files", Files may range in size from 768-bytes up to the 
limit of the magnetic tape storage capability. Files may be 
either ASCII or binary in format. All files in the 
Operating System are stored in Binary format to take 
advantage of a peripheral “firmware package called the 
binary loader. inary data transfers between the processor 
compiler and magnetic tapes are faster than when stored in 
ASCIY. This is because the processor uses binary format for 
internal processing and the conversion from ASCII to binary 
format is eliminated when data is stored binary. 
Additionally, a stored binary program can be transfered to 
the processor memory from a specified peripheral device 
without disturbing variables and associated values 
previously defined. This feature is not available to ASCII 
programs. 

The GS permanent ROM has been programmed to respond to 
BASIC (Beginner's All-Purpos2 Symbolic Instruction Code), a 
high-level programming language first developed at Dartmouth 


College. The Operating System has been correspondingly 


97 





written in this language. However, BASIC as used by the GS 
has been extended to fit the language to the specialized 
capabilities of Graphic System. These extensions are in the 
areas of graphics, file system access, unified handling of 
input/output operations, Matrices, charecter Steang 
Manipulation, high-level language interrupt handling, and 
Operating system facilities. A key difference between 
extended GS BASIC and most other BASIC languages is that 
most keywords and their parameters can be evaluated 
independantly of program control. The result is a Lrich, 
versatile, yet simple to use programming tool. 

Three peripheral ROM firmware packages have been added 
to the TDRL system to tailor it to specific needs of the 
laboratory. may ace the Signal Processing ROM Pack No. 1, 


Signal Processing ROM Pack No. 2 (FFT), and the EDITOR. A 


description of each is provided in Appendix A. 


C. OVERVIEW OP THE TDRL OPERATING SYSTEM 

The TDRL Operating System has been specifically designed 
+o be higly operator-processor interactive. This provides 
the greatest degree of flexibility in the acquisition, 
processing and evaluation of electromagnetic time domain 


transients for target imaging. 
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The total Operating System occupies 83818-bytes of mass 
storage on TDRL Library Tape 1 in three files. It requires 
an additional 201628-bytes of storage in 62-files on TDRL 
Library Tape 2 (Data Storage). The Operating System is 
logically divided into three-individual programs: 
1. INPOT--for the acquisition and averaging of signals 
of interest. 

2. MATH--performs the significant mathmatical processing 
of the acquired signals for inverse scattering. 

3. GRAPH--provides a visual display to the operator of 


the results of the processing. 
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Mass storage requirements of the Operating System 
programs are given in TABLE III. The storage requirements 
listed in TABLE III are for the programs only and do not 
include requirements for the necessary definition of arrays, 


Strings or numeric constants and variables. 
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All programs are designed to be as nearly independant of 
each other as possible, after an initial run, to allow for 
naxinum FPexibility and Weeiety in data usage and 
2valuation. The most dependant program is MATH which cannot 
be used without inputs from either INPUT or GRAPH. The most 
independant program 1S GRAPH which is capable of being run 
with inputs from mass storage only. 

A general flow diagram of the operating system is given 
in Figure 4.2. Specific named subroutines have not been 
shown. The next section will describe in detail the 
individual programs of the Operating System and their 


specific requirements and processing methods. 


D. OPERATING SYSTEM PROGRAM DESCRIPTION 

This section will describe in detail the individual 
programs and their subroutines that make-up the Operating 
System for the TDRL. A general description of the purpose 
of each program will be provided first. Sources of inputs 
and outputs will be noted. A general flow diagram for the 
specific program will be provided Showing how all 
Subroutines interact. This general overview will then be 
followed by a more detailed description of each subroutine 


in the progran. The purpose of the subroutine will be 
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noted. Sources of input and output will be defined. A 
discussion of applicable algorithms and the imeLfewti (ene 
information will be also provided. Zxamples of visual 
displays applicable, potential problems and correction 
methods, will be discussed when warranted. 

As noted, the following commentary will make frequent 
referance to program generated messages. Examples of these 
messages will not be given in the pvody of this text, but can 
be found in Appendix F, a general run of the Operating 
Systen. 

1. Program INPUT 


INPUT provides the means to initialize and drive the 


total Operating System. This sub-program is stored on 
Library Tape 1, File 2. Memory space required for its 
implementation in the processor CPU is 24552-bytes. Nase n 


all arrays and constants defined, a total of 42469-bytes of 
memory space may be required. A listing of INPUT can be 
revi2wed in Appendix C. 

INPUT is the front-end program for the Operating 
System and as such, will generally be the first program to 
be initiated for TDRL measurements. Program start-up is 
begun by placing Tape 1 into the 4052 built-in magnetic tape 


drive, and placing Data Storage Tape 2 in the 4924 Remote 
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Tape Drive. INPUT is then normally placed into the CPU's 
working memory under operator control by the issue of the 
following commands from the keyboard: 
FIND 2 
CALL "BOLD" 

Program execution is next begun by typing the command "RUN" 
from the keyboard. 

If this is the first attempt to run the Operating 
Systen in conjunction with the Digest a ll Processing 
Oscilloscove (it is assumed the DPO is "ON" whenever INPUT 
is employed), the following message may appear on the GS CRT 
BoOuLOWing the issuance of “FIND 2": 

“No SRQ on anit in immediate line -- message number 43." 
Iz the program is being run from any line prior to line 140 
(the usual case), this message may be ignored and the 
program will automatically correct the e2rror message and 
continue processing. If the program is being run subsequent 
to line 140 (the case when returning to INPUT from some 
Other internal file), the following command mus* be entered 
from the keyboard: 

POLL N,M31 
Following this, the regular sequence of conmands may then be 


issued to start the systen. 
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As noted, INPUT aay alternately be loaded and 
mmeeclLated under program control from GRAPH (File 4&4). 
Execution will begin following line 180 in this case. This 
return to INPUT occurs in normal Operating System (OS) 
Operation after an initial run has been completed and 
information is then desired for an entirely new target. 

The program specificaily provides for the means to 
input information concerning the transient electromagnetic 
response of a target. It acquires, averages and stores all 
pertinent waveforms. INPUT also allows for the easy storage 
of specific antenna and target dimension parameters for 
comparison when the target is known exactly. This data is 
permanently maintained in a library on Tape 2 in Files 5-49. 
INPUT further initializes most of the parameters used later 
in Program GRAPH for the visual display of the processed 
results. 

Data inputs are from the keyboard or internal 
memory. Outputs are to internal memory, mass storage and 
Bite 3 (MATH). INPUT is highly operator-processor 
interactive. Numerous directives and visual displays guide 
the user. 

the structure of the progran consists of 


11-subroutines. One (DRIVER) is of general application to 
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the progran. Five (ACAVG, INWAV, DIWAV, AUGWAV, SCATWAY) 
are specifically related to input, acquisition, averaging 
and storage of relevant information concerning the 
time-domain transient response of the target in question. 
The remaining five subroutines (ANTLOC, IGTLOC, NEWPAR, 
ANTIN, TIGTIN) provide special services for ae operating 
Systen. A general flow diagram for input is provided in 
Figure 4.3 A description of each subroutine follows. 
ae DRIVER 

Mises the Main driving routine for INPUT. ie 
consists of 129 steps (40% of the program) from line 100 
through line 1260. It is initiated directly under operator 
Sen~rolat tine 100, or under program control from GRAPH in 
line 180. Data is input from the keyboard or from the mass 
storage files located on Tape 2. DRIVER will pass program 
control directly to MATH when directed. 

Lines 100-300 initialize various variables and 
dimension various arrays. Line 150 clears the General 
Purpose Interface Bus (GPIB) of the DPO service request. 
Plag ?7 (indicates whether the acquired data is destined for 
use in Time Domain or Prony Methods) is set to 0O (Tine 
Domain) in line 160. P8, the run number, is initialized to 


zero in line 170 and will be incremented by one in line 210 
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Beee every ran (including the first). Various string 
constants are initialized in line 240-250 and 280-290. Plag 
P9 (indicates the need for initializing input parameters) is 
set to 1 (input parameters reguired). Lines 320-530 provide 
a menu of subroutines and functions which are available for 
Operator selection. The selection of a process from the 
menu is input from the keyboard into system variable PQ. An 
example of the menu is given in Appendix F. 

Lines 550-570 direct the program flow to the 
appropriate portion of DRIVER. if the operator has selected 
one of the first four subroutines, and this is an initial 
run of the 0S, then the time for a full screen sweep of the 
DPO must be inputted to the program from the keyboard. The 
sweep rate is put in in seconds in T1 in line 640 and must 
be greater than or equal to the time reguired for the the 
DPO beam to sweep completely across the DPD screen. 

Program control is now passed to line 850-1070 
where the following cun descriptive parameters are input: 

1. Date 

2. The distance the target is from the transmitting 
antenna 

3. The antenna number selected 


4. The target number selected 


NO, 








5. Explanatory remarks 
These are all keyboard entries. An exampl2 of the messages 
that query entry of the above parameters is given in 
Appendix F. 

ier Gane as not the first for the Operating 
System, a different series of inputs are queried. En eas 
case, program control is passed to lines 670-840 which 
generate messages on the GS asking the operator if there is 
geeecnange in the original location of the target, if a new 
antenna 1s to be used, or if adefined target is to be 
changed. The operator inputs the alpha character "Y" or "Nt 
por "yes" or "no" fromthe Keyboard to reflect the current 
Situation. Appendix F displays these messages. 

Lines 1000-1040 generate the input of the 
antenna parameters. These parameters ar2 drawn from files 
stored on Data Tape 2. Each antenna file, located in one of 
the files numbered 5-24, is 2 63-element array, returned in 
AS. In elements 1 through 60 are located: 

lee emecUN number $£Or this particular antenna-target 
combination 
2. the antenna length compensating parameter for the 


antenna-target combinaticn 
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Bae che antenna amplitude compensator OG the 


antenna-target combination 


AS can pe seen, antenna parameters are actually three 
elements within the larger array, each three elements 
identifying a unigue antenna-target combination. The 


parameters in elements 61-63 are system antenna parameters 


used for unknown targets. They are sinply the sum of all 
run, antenna length compensator updates and antenna 
amplitude compensator updates. The average of the last two 


elements is applied to. signals acquired if the target is not 
defined from the target library. This is shown in lines 
1024-1025. 

The antenna amplitude compensator is a very 
sensitive function of the distance the target is from the 
transmitter due to EM propagation attenuation. In line 
1025, for unknown targets, and in line 4152 for defined 
targets, these losses are taken in consideration by finding 
the square of the ratio of the targets true location on the 
plane with referance to the transmitting antenna and the 
control distance (1.27m). 

After the definition of variables is complete, 
DRIVER passes control to the appropriate function in line 


1090. Line 1100 returns the display to the Menu for further 
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function selection or to line 1170 wher2 the DRIVER passes 
Mencsel to MATH and terminates any further INPUT processing. 

In summary, INPUT DRIVER reguires input of 
selected program routines and initial parameters by the 
operator from the keyboard, and requires the 0S to determine 


whether the current run is the initial for the day ora 


Subsequent one. The parameters are saved for later graph 
requiremts. Funtion selection drives the appropriate 
Subroutines and routines. A flow diagram for GRAPH DRIVER 


is provided in Figure 4.4. 
be. Subroutine ANTLOC/TGTLOC 
ANTLOC and MTGTLOC are Similiar subroutines in 
that they locate the appropriate antenna and target storage 
arrays filed on Tape 2. Antanna parameters are located in 
Piles 5-24. Target parameters are located in files 26-45. 
Lines 4000-4060 contain ANTLOC. TGTLOC is found in lines 
4100-4160. HOreoc has the Madditional Eunction of defining 
the system antenna parameters for the current run, as shown 
in lines 4149-4153. Ther? are no operator-procecessor 
‘interactions during execution of these subroutines. 
c. Subroutine DIWAV/INWAV/AUGWAV 
DIWAV, INWAV and AUGWAV are sinilar subroutines. 


They drive ACAVG to input the direct waveform, the incident 
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wavefora, and the augmented waveforn. They allow the 
selection by the Operator of the number of signals to be 
acquired for each waveform by the DPO (total composite 
waveforms from the sampler) and averaged by the OS. Each of 
these subroutines is accompanied by a nessage on CRT 
describing the procedure to be employed in acquiring the 
Bacticular Signal. These messages are reproduced in 
Appendix F. DIWAV is located at lines 3000-3170. It saves 
data in File 50 of Tape 2. INWAV is located at ines 
3200-3380. It saves data in File 51. AUGYAY is located at 
lines 3400-3590. It saves data in File 52. 
d. Subroutine ACAVG 

ACAVG acquires and averages the appropriate 
Waveform interactively requested by the operator from the 
processor keyboard. It also removes any dc signal that 
might be included in the acguisition. This feature allows 
the operator to position the desired Signal anywhere2 on the 
DPO screen and still maintain a true zerd axis in later 
Seaphic displays. Note, the signal waveform must be 


completely within the DPO screen area. ACAVG is 74-steps 


long (23% of the progran) and is located in statements 
2000-2730. It is driven by one of the appropriate input 
waveform subroutines previously described. Fo Gecvras 
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results to that subroutine for later storage. It also 
displays system information for the current signal being 
acquired to the DPO screen. 

Tznes 2090-2130 initialize ACAVG by requesting 
the operator to input the number of waveforn signals to be 
averaged in RZ This value is then returned for 
verification. In lines 2140-2540, the OS writes the current 
acquisiton number of the waveform to the DPO screen, 
acquires and stores the Signal into DPO memory and then 
transfers the signal to the 4052 for further processing. 

Lines 2170-2370 perform the count function. 
Depending on the number of waveforms desired, the count will 
be for every signal (N2 less than 21), for every five 
Signals (N2 less than 101) or for every ten signals (N2 
greater than or equal to 101). This 1s done to reduce the 
time required for acguisition 9f all signals. 

Lines 2410-2430 instruct the DPO to store and 
hold the waveform in DPO memory. Line 2420 stops processing 
to allow for the DPO to fully acquire the waveform ina 
Single sweep. T1 is the length of the wait. 

In lines 2440-2050, the GS instructs the DPO to 
transfer the digitized signal from DPO memory to the 


microprocessor temporary storage array X0. Lines 2470-2500 
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then search the first eight values of the array for any 
obviously out of range values input by poor DPO blanking on 
the signal front end. If such a value is found, the array 
element is set to discrete value 511, corresponding toa 
scaled value of 0. The waveform is next added to the sum of 
the previous waveforms for the current acquisition in line 
2550. The resultant matrix is the time sampled values of 


moe signal in BO: 
BOP ee XO” (4.1) 


If further acguistions ar@ necessary, th2 entire process 
described above is repeated. LienoO Luccthss acquisitons are 
to be done on the current signal, the ensemble average tine 


Sampling matrix xX0 is determined in line 2550: 


(4.2) 
XO = BO/N2 


Any de offset is eliminated in line 2560. 

It is interesting to note the tine required to 
acquire and average signals using the solid state UHF pulse 
generator. For ten averages of a signal, 42 saconds of 
elapsed time is needed, or about 4.2 seconds per acquistion. 
Of this, about 1 second is needed to write information to 


the DPO screen per acquisition. For 21 averages, about 65 





seconds of acquisition time is required, or about 3.1 
seconds per acquisition. Thasails about a 35% reduction in 
acquisition time overall per signal. The improved time per 
acquisition is due to the reduced frequency with which data 
in the form of the wavecount is written to the DPO screen. 
For a relatively small increased investment of time, a 
tremendous advantage is realized in averaging signals and 
reducing polluting noise. A further decrease in acquisition 
time per signal is realized when 101 or nore averages are 
specified. However, as will be discussed in Chapter V, lack 
of total signal stability may produce errorneous results for 
such a large number of acquisitions due to the length of 
time required. 

In lines 2570-2820, the waveform scale factors 
are determined and transferred to the processor. The values 
are returned in M$ (vertical magnitude scale factor) and T$ 
(horizontal time scale factor). THe discrete,  digitzzed 
verticle values are appropriately scaled to the correct real 
tine values in line 2700. The time between a single 
discrete sample on the horizontal axis (sampling rate) is 
returned in Z1 at line 2820. 

In summary, ACAVG does the physical acquisition, 


averaging and scaling of each signal. It returns the 








current waveform count to the DPO CRT for display, and the 
appropriately scaled waveform in XO == ene! Ginseng 
Subroutine for storage in the appropriate file. fe fuGeuer 
returns the appropriate horizontal and vertical axis labling 
scale factors in T and M respectively. The flow diagran 
for ACAVG is provided in Figure 4.5 

The basis for this subroutine was derived by 
Capt. Hammond. [Ref. 19] His program in its entirety has 
been extensively modified to fit the current requirements of 
ene [DRL OS. 

@. Subroutine NEWPAR/ANTIN/TGTIN 

NEWPAR, ANTIN and TGTIN are interactive 
Subroutines which allow an operator +o input new parameters 
for defined antennas and targets into the mass storage 
library on Data Storage Tape 2. 

NEWPAR is the driving routine for ANTIN and 
[GTIN. It receives its input from the DRIVER and selects 
either ANTIN or TGTIN. It returns the new parameters to the 
appropriate file on the storage tape. NEWPAR then determines 
if further inputs are reguired. fe evOt, OLOgGEamM GCOntEOl eas 
returned to DRIVER. 

ANTIN defines the values of the parameters of 


the new antenna. P= Oiosts stne Operator to input the 
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antenna number. An array of 63 elements is then generated, 
each third element being filled with an initial antenna run 
humber equal to 1. All other elements are set to Zero. Et 
the operator specifies an antenna value greater than 20, the 
following message is returnei. 

"STORAGE AREA LIMITED TO 20 ANTENNAS. WANT TO CONTINUE?" 
This is necessary as current space limitations limit the 
humber of antennas to be defined. 

TGTIN defines the values of the parameters of 
new targets. It requests the operator to input the new 
target number as well asthe length and radius of the 
earget . The carraige return must be pressed after each 
en cLy. 

NEWPAR, ANTIN and TGTIN are located in lines 
4200-4800. They comorise 11% of the progran. 

fF. Subroutine SCATWAV 

SCATWAVY determines the target backscattered 
waveform based on the differance of the augmented and 
incident waveforms. It returns the backscattered waveform 
Momcugsagemat tile 53 On Tape 2. It also sets flag P? to the 
appropriate value required for either Prony processing 
(P7=1) of Time Domain processing (p7=0). Program control 


1s then returned to INPUT DRIVER. 
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MATH performs most sf the significant mathmatical 
computations for the 0S. its cerrect performance is central 
in the production of 2ccurately processed signals. MATH is 
stored on Tape 1, File 3. It reguires 14544-bytes of CPU 
memory space. With all arrays and constants defined, a 
maximum total of 49108-bytes of memory are required. 

This sub-program is initiated under program control 
only from INPUT and GRAPH. There is no interaction between 
operator and system during MATH's execution. 

Inputs to MATH are from INPUT, GRAPH or Tape 2. 
Values are returned to mass strorage or GRAPH. 

MATH performs the following functions for the OS: 

1. Past Fourier transforms the Direct and Backscattered 
waveforms. 
2. Determines the frequency domain impulse response. 
3. Performs inverse Fourier transform to find the tine 
domain inpulse response. 
4. Determines the step and ramp response for the target. 
5. Determines the physical optics shape of the target. 
6. Optimelly tiiters the fraquency domain impulse 
response, yielding a smoothed time domain impulse 
response, ramp response and physical optics shape of 


the target under noisy conditions. 
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7. Compensates the ramp response and optics shape for EM 
propagation attenuation. 

8. Updates optimal antenna parameters for the antenna 
being used, based on the targets ramp response, 
automatically. 

MATH is composed of four subroutines and DRIVER. A 
description of each follows. A general flow diagram for MATH 
is given in Pigure 4.6. 

as MATH DRIVER 

The driver routine for MATH is located between 
lines 100 £0 1210 and comprises 59% of the total progran. 


In addition to initializing the program and driving the 


subroutines, DRIVER performs most S59 the -Ssidgnmayeicant 
mathematical computations in straight line data flow 
fashion. 


MATH DRIVER is initiated at line 100 from INPOT 
Be lire’ 130 £Erom GRAPH. Parametric data is input directly 
from INPUT. It is recovered from mass storage in line 160 
Pe GRAPH initiates Ward. Lines 100-221 serve to initialize 
the local program parameters. 

Line 250 inputs into x0 from mass storage file 
53 the signal (backscattered waveform) which represents the 


convolved transient response of the transmitting antenna, 
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the receiving antenna, andthe target. This data is then 
used to drive subroutine FFT. Oe se returicd feom rr?! Jn 
and P as the frequency domain magnitude and phase and saved 
in arrays E and 53 respectively. 

A similar procedure is performed in lines 
290-340 for the convolved transient response of the 
transmitting and receiving antennas without target on the 
ground plane. The magnitude and phase are returned in 
arrays F and F3 respectively. 

In lines 360-420, the frequency domain impulse 
response of the target and noise optimization are performed. 
loni tude of the impulse response is returned in £0, with 
phase returned in FO. The algorithm for determining the 
Lmupulse response and performing optimization were developed 
in chapter III and will not be repeated hers. 

Lines 440-540 perform the inverse Fourier 
transform of the frequency domain impulse response of the 
target to find its time domain equivalent. The frequency 
domain data entering at line 440 is in polar form. This 
form, is most suitable for the previously described 
mathematical manipulations. It is not a suitable fora for 
the inverse Fourier transform nethod employed by the 4052 


system. In order to perform this operation, the magnitude 
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and phase arrays must first be converted into two source 
arrays containing real and imaginary data. Next, tne EFXrst 
element of the real array must contain the signal dc tern, 
while the last real element must be the Nyquist frequency 
value. As both these elements represent pure real numbers, 
ene first and last elements of the imaginary source array 
muse be correspondingly zero. 

Lines 440-470 establish a FOR/NEXT loop that 
returns real terms to M and imaginary terms to P. Real 
terms are established by multiplying the polar magnitude, 


EO, by the cosine of the corresponding phase FO: 


M(j) = E0(j) x (cos (FO(j)) (4.3) 


Imaginary terms ar2 established by multiplying the polar 


Magnitude by the sine of the corresponding phase tern: 
P(j) = E0(j) x sin(FO(j)) (4.4) 


INPUT had previously removed ali de levels prior 
to MATH, so M(1) and 4(257), the dc and Nyguist terms, are 
Set aqual to zero in lines 480-490. Lines 500-510 set the 
corresponding imaginary terms to zero. Note that both of 
these arrays must be 257 data points long, as the output 


feed yeowoet be 512 data points. 





Data is now correctly formatted in the two 
source arrays so that the command "“INLZSAV" can interleave 
the real and imaginary data into a single destination array, 
10, whose format is acceptable to the GS "IFT" command. 

Line 540 calls the system inverse Fourier 
transform command "IFT" which transforms the Spectral data 
into waveform (time-domain) data. The spectral data 
originally in I0 is overwritten by the new time domain data 
aod returned in I0. TO is now the time domain impulse 
response. 

The inverse Fourier transform can be expressed 


mathematically by the following summation: 


es j 2mnk/N 
a a7 Ny re se corer =08 eel cans) 
In the above equation, N refers to length of the array 


argument, andn iis an index used in generating the various 
elements of the output array. X (k) is the k-th conplex 
Fourier coefficient, and X(n) is the n + 1-th element of the 
real data output. (Ref. 20} The 4052 uses an extremely 
Paro & method for Ea computations, the Sande-Tukey 


decimalization-in-frequency algorithn. 
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Line 620-720 save the time domain impulse 
response (file 54), and frequency domain direct waveforn 
fete 55), the backscattered waveform (file 56) and the 
impulse response (file 57) in mass storage. 

Next, the time domain step response of the 
impulse is returned to X9 in line 740 by using the 4052 
"INT" routine. This is an integration calculated using the 
trapezoidal rule for approximating the definite integral as 


follows: 


B (t) = B(t-1) + .5 * (A(#-1) + A(t)) (4.6) 


SOL eee = -2oedigeciete pil 


where; A = the source array 
B = the destination array (result) 
N = the number of elements in the array 


Any de value in the result is removed through the actions of 
lines 750-760. 

yeas , the time domain ramp response is 
determined by integrating the step response in X0. The ramp 
ls returned in RO. The MATH DRIVER outputs the time domain 
ramp and its minimum value in R2 to Subroutine OPTICS which 
returns the physical optics shape (shape) and ramp response 


compensated for electromagnetic propagation attenuation. 
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Note the shape does not aS yet have the correct amplitude 
values in its array elements. 

BEtnewgin 1S the Elest for the target (prior to 
optimization) and the target is one defined in mass storage, 
MATH DRIVER inputs the last computed antenna parameters for 
the compensator for shape length (S5) and radius (S6), as 
well as the actual defined lenghth (C6) and radius (C5) to 
Subroutine ANTUP. SS and S6 are then automatically updated. 
The physical optics array is then multiplied by the 
compensator value S5 so that its elements reflect the 
estimated dimensions of the targets in inches. This is 
based on the experimentally observed fact that a particular 
size and configuration of stationary target will return some 
backscattered field relative to its size and shape. That 
this value is not exactly the same from one run to the next, 
but only approximately so, is due to certain stastically 
determinable effects produced by noise, small differences in 
Seaentation, etc. A large number of runs should produce an 
averaging value which takes into account these factors. 138 
the target is unknown (the general case), or the current run 
is for noise reduction, previously computed values of the 
antenna parameters are recalled from mass storage. Por the 


unknown target these compensators are the averages for all 
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meeoets 2<Or all runs fora particular antenna. ia Lace, 
OPTICS is not part of the data processing flow, the 
compensators having been previously deternined in INPUT at 
lines 1024-1025. 

Following a known target run, the moininum value 
@eepe 2S Seturned to B8. If the run is not for optimization 
(L1=1), then B36, an amplitude compensation cptimization 
factor is set equal to B8 and the actual length of the shape 
me retained in Q9. For graphing purposes, it is desirable 
that B6 and Q9 not vary from optimization eb, awe, 
Seeregization run. MATH DRIVER will input BO, Q9, and B6 to 
Subroutine OPTIM, which compensates the amplitude and length 
of the physical optics shape for the low pass filtering done 
to reduce noise. An arbitrary target is treated in exactly 
the same manner. (Note, this effects the shape ONLY and not 
the time domain ramp response during optimization.) 

The above system has worked remarkably well in 


deducing the actual physical parameters of both known and 


unknown targets. *®Or known targets, accuracy is invariably 
foe OT less. eoue unknown targets, this ‘accuraey 
deteriorates to about 15% or less. If the data is treated 


in the raw form, with no compensation applied, accuracy is 


metieeca tO 25% O05 worse, up to 50% for the unknown target. 


i 





This can be noted on the graphic output of the shape where 
the compensated, uncompensated and known (when given) values 
See Length abe compared run to run. 

Lines 1900-1140 save the updated antenna 
parameters in the appropriates files, the shape array in file 
BIO» the time domain ramp response in file S59, and 
Significant system parameters in file 60, all on mass data 
storage Tape 2. 

Lines 1160-1190 is MATH DRIVER ending routine, 
deleting variables and passing control of the program to 
GRAPH (file 4) on Tape 1. 

In summary, MATH DRIVER is the principal vehicle 
for OS processing for inverse scattering. A flow diagram is 
given in Figure 4.7. 

b. Subroutine FFT 

FFT is the fast Fourier transform subroutine 
used tc convert the time domain Direct and Backscattered 
waveforms to spectral data. This subroutine is located in 
lines 2000-2140 and requires 7% of the total program space. 

XO contains the real-valued signal data on which 
the Pourier transform is to be perforned. It is input from 
Mark DRIVER. Line 2060 calls the GS command ‘rhe sand 


performs a fast calculation (about 4% seconds) of the DFT 
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(discrete Fourier Transform), expressed nathematically by 


the summation: [{Ref. 21] 


1 sodeees 
X4(k) = 7 Maver” | cor kao 1... N/2 (4.7) 


where: N = array length 


k = index used in generating various 
Fourler coefficients 


k-th Fourier coefficient 


Ky (k) 


X (n) n+ 1I-th elemant of the real data 


input 

After the command has executed, the original 
real data will be overwritten by the results of the FFT 
computation. The output data is in complex rectangular fora 
With the first element of the array X0 containing a real 
Rumber representing the dc value of the signal, the second 
element containing a real number representing the value of 
the discrete Fourier transform at the Nyquist frequency, and 
the remaining array elements containing alternately real and 
imaginary Fourier coefficients. 

This arrangement 2S “nee suitable for the 
division that later takes place. In line 2070, the 
rectangular array X0 is converted to polar form by calling 


the GS command "POLAR". Magnitud? components are returned 
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in array MH and phase components are returned in array P. At 
mars pOlnt, SUDEOUtIne FFT returns control to MATH DRIVER 
Meter ne abbays XO, M and P for Eurther processing. 

Cc. Subroutine OPTIC 

Subroutine OPTIC manipulates the ramp response 
of the target to recover the physical optics shape, but it 
does so without effecting in anyway the raw data. It later 
compensates both th? ramp response and the shape for EM 
propagation attenuation experienced along the body of the 
target. Subroutine OPTIC resides in lines 2200-2700 and 
comprises 25% of Program MATH. 

The ramp response RO andthe location of its 
Minimum value within the array are input into OPTIC from 
MATH DRIVER. The mean value of the response is calculated 
in lines 2250-2260 and returned in R3. This value is used 
later in determining the endpoints of the shape array. 

Lines 2290-2500 is a FOR/NEXT loop that 
seperates out the targets shape from the ramp array data. 
The search for valid data begins at the point in the array 
equal to the location cf the mninimum value. It is assumed 
that the target will always return a response larger than 
any that might be due to range noise or other excitation. 


For the size of the targets being used, and the the 


131 





configuration of the present source, this has always proven 
+o be an accurate assumption. The FOR/NEXT loop searchs 
first for a tzrailing-edge andpoint. This is done in 
ascending time steps in search of the point where the array 
elements become zero or poSitive, representing a zero-axis 
SEossing. For every value of RO less than zero, a 
corresponding value is input in BO, the shape array. ila 
BO(k) is a positive value, representing a zero-crossing, the 
FOR/NEXT loop is terminated. If BO(k) is negative, OPTIC 
then checks the location of the point. it Lsaimport da ltmro 
do this as noise and instability of input signals from an 
acquisition-to-acquisition run may erroneously shift-zero 
crossings or preclude zero-crossings altogether. In these 
cases, determination of the shape must be somewhat nodified 
as described below. 

meeetan> 2350, two situations are checked. 
First, is the preceding array element of RO more negative 
than the present element and is the present element greater 
that 90% of the cranp cesponse average value. ise 
affirmative, the next array element is fetched. A positive 
result of this check indicates that datum is on the skirt of 
the ramp response, and that the data magnitudes are 


progressing towards a zZero-axis crossing normally. If the 
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Check returns as negative, that is, the present array value 


is more negative than the pr 


Ww 


ceding slement, than an anomaly 
has occurred either due to noise, or time-axis shifting. In 
this case a special endpoint routine must be employed. 

The second check is to determine if the array 
element is oon the main Signal portion of the response and 
not the skirts. If the element amplitude is greater that 
90% of the mean value, than it is assumed that the element 
1s on on the skirt irregardless of the reiationship between 
the amplitude of the present and preceding elements. The 
next datum is fetched. It can be seen that this arrangement 
ls subject to some error if a nolsy el2ment is encountered 
just at the point where the ramp response passes the 90-th 
percentile ~~ the signals mean value. This error is 
considered relatively minor. 

The 90-th percentile was chosen foliowing 
numerous observations that noise or time-shift in the main 
ramp response normally do not produce fluctuations that 
cross this boundary. Error in the shape representation will 
be introduced if it does, in fact, cross this levei. The 
result will be a fore-shortened drawing of the target. 

LE ».the need for an endpoint routin? is 


indicated, processing goes to lines 2360-2470. In line 
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2360-2380, the distance the last valid data element 
(RO (x-1)) is from the zero axis is determined. In line 
2350, the distance between the previous two valid data 


elements is found and returned in W(3): 
W(3) = RO(k-2) - RO (k-1) (4. 3) 


The above is for the ascending time series run of the 
FOR/NEXT loop (trailing edge). For the descending tine 
series, the values in the parenthesis have opposite sign. 
This difference represents an approximation of the slope of 
the skirt. A problem with this approximation is that it 
often occurs where the slope is not as great as the average. 
The result 1s a somewhat elongated figure. + may be 
possible to defeat this problem by taking more than two 
points for the average. 

In lines 2370-2380, a check is made to see if 
the last valid data point is within 10-times or less of the 
difference determined above. Le Met de Sie the actual 
difference is returned in Q71 and used +9 drive the values 
Beaced in BO to a z@ro-crossing. T= the difference is less 
that 10% of the amplitude needed to effect a zero-crossing, 
then one-tenth of the actual distance is returned in Q1. [In 


this manner, the shape rapidly goes to zero and an endpoint 
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Ss determined close to the estimated real endpoint hada 
Natural zero-crossing occured. The endpoint is located with 
a FOR/NEXT loop between lines 2420-2470. This FOR/NEXT loop 
Simply iterates the values of BO by Q1 until a zero-crossing 
is effected. Whenever BO becomes positive, all loop 
processing 1S terminated and the endpoint is returned in 
W(1) £or the trailing skirt and W(2) for the leading skirt. 

If the above process was for the leading edge of 
the ramp response, the loop is reentered and the trailing 
edge check of the data is performed. The method employed is 
identical with the above process for the leading edge. Note 
meet up to this point, the ramp response in RO has been 
totally unaffected by any processing in OPTICS. It has been 
used simpiy as source of data. 

Once the target tesponse has been located in RO 
and entered in BO, then both arrays can b2 compensated for 
propagation losses experienced. This is done in lines 
2570-2640. The value of the nominal distance the target is 
from the transmitting antenna is returned in Z6. A FOR/NEXT 
loop in lines 2590-2640 corrects the arrays itor the 
estimated propagation losses. Line 2600 computes the 
distance along the body the wave propagates for each sample 


Beant in RO: 
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R = (Q - (1) * 21 * (10*8) /2 (41-9) 


where: (9 - W({1)) * Z1 = twice the time of travel 


of the propagated wave 


Pa line 2610, the propagation loss compensation factor is 


returned using: 


ee Zone e+R) (4.10) 
ae 


Where; Z = range from transmitting antenna to 
wasgec face 

It is also assumed that the distance between transmitting 
and receiving antenna is Z. Z is in fact the measured value 
of distance read from the "TIME/DISTANCE" dial on the DPO. 
Bach value of RO and BO is then multiplied by this 
Gempensation factor in turn. Processing continues once all 
values for the computed length have been approoriately 
amended. 

Line 2260 is a FOR/NEXT loop that transfers the 
values ln the array BO to the beginning of the array to 
facilitate graphing. SaeCewrnls LS accomplished, progeanl 
control is returned to MATH DRIVER with the shape unscaled 


for actual dimensions of the target (values are equal in 
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amplitude to the compensated voltage magnitudes of the ramp 
response). The ramp tesonss 1S returned uneffected by any 
of the above processing except for propagation attenuation 
compensations. Figure 4.8 provides a flow diagram for 
mmoroutine OPTICS. 
d. Subroutine ANTOP 

This surboutine updates the specific antenna 
length and radius compensation values if the target is one 
defined from storage on Tape 2. If the target is arbitrary 
in the general case, no updating occurs. ANTUP is located 
between statements 3000-3080 and comprises 4% of the 
program. 

The unscaled shape array BO is input to ANTUP 
me Om MATH DRIVER. In line 3060, the amplitude compensating 


and scaling factor is computed using: 
DO= (A + S*B)/(C*B) 


where: S6 = scaling/compensating amplitude factor 


‘2 
l 


= radius of target 


8 = minimum magnitude of ovhysical optics 
array 
RK = the number of previous times this 


antenna/target combination was used 


S = sum of previous compensation values 
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The result is sinply the averaged value of 
Sompensation/scaling for allruns using this varticular 
antenna and target combination. 

Line 3070 performs a similar computation for S5, 
updating the compensation/scaling factor for the target 
imength. The formula is identical except that the defined 
target length replaces the target radius and the observed or 
measured length from the plane replaces the mininunm 
Magnitude of the shape array. The concept is that slight 
differences in target returns due to variations in target 
location and orientation, noise, etc., will be averaged out 
over the long term for any particular antenna. When the 
target is arbitrary and there is noa priori information 
then ANTUP is bypassed and the average of all the 
compensation/scaling factors for all targets for the 
particular antenna is used. This requires that any new 
antenna be "seeded" with compensation/scaling factors from 
at least one run uSing known targets prior to attempting to 
run an arbitrary target on the antenna. What results is a 
"poor man's" compensator for the currents on the body of the 
target based on empiracal observations of many target 
responses. It is not as accurate as the full solution of 


the tine domain integral equations would be, nor is the 
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exact target shape returned. What is obtained is an 
excellent estimate of the target, whether known of unknown, 
Desed On an historicai compilation of data for the TDRL 
system. This concept can be generalized to any system using 
components which are nominally fixed, such as large scale 
bistatic radars. The shape itself, az least for simple 
axisymmetric targets, can probably be returned with the 
development of the proper software routines. Perhaps the 
rate at which the ramp response slopes vary may be analyzed 
Ber this purpose. What is obtained is a method which 
employes a small scal2 processor requiring relatively simple 
processing to obtain accurate time domain estimates based on 
the ramp response of the target almost in real time, and 
capable of being generalized to almost any system employing 
the proper source. 

The compensation/scaling factors are returned to 
MATH DRIVER along with program control. BO is unaffected as 
We= by SS or S6. 

e. Subroutine OPTI 

OPTIM compensates the shape array for the 
effects of low-pass filtering of the target responses using 
Riad's Optimization technique. Optimization is done to 


remove high-frequency noise terns and, for the Snape, to 
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improve the overall shape of the drawing. ie cea lacy, 
without compensation, the drawing is less representative of 
the true targets shape when optimized. Low pass filtering 
Will lengthen the overall shape and reduce the radius, the 
amount being relative to the value of lambda, the 
Optimization variable. OPTIM allows for this and introduces 
factors which will pin the optimized shape to the original 
form while allowing the benefits (reduced noise ripple on 
the body) of low-pass filter noise reduction to be shown. 
OPTIM is located in lines 3100-3180 and represents 4% of the 
progran. 

Bo, 86, (original minimum of the unoptimized 
Shape), 88 (the present minimum of the shape) and Q9 (the 
Oliginal unoptimized target length) are input from MATH 
PReVER. The amplitude OPTIM compensation is found in line 
Bamyo by dividing B86 by B8. It is returned in B7. The new 
compensated target length factor is returned in Q7 and is 
found by dividing the present value of the length py the 
Original value. Finally the whole physical optics shape 
array is scaled by the new amplitude factor. The new values 
are returned to MATH DRIVER, along with the original factors 
in B6 and Q9. menact, arter the initial unoptzimized run, 


B6 and Q93 are constants until a new target is acquired. 


141 





3. Program GRAPH 

GRAPH provides a central means to visually dispiay 
the results of processing performed in MATH . VieneOnRTneCet 
evaluation and analysis of the acquired processed data is a 
merect function of the effective reproduction of the 
resultant arrays by GRAPH into a form nost useful and 
understandable to the operator. GRAPH is located on Tape 1, 
Pile 4, and requires 44712-bytes of memory space. 

Program GRAPH is normally initiated under program 
control from Program MATH, in line 100. It may be initiated 
directly by the operator if the results of a previous run 
have been retained in mass storage. The commands issued by 
the operator are; 

FIND & 
CALL BOLD 
RON 
Meeeputs are to the CRT or a return of control to one of the 
other vrograms that make-up the TDRL OS. No results are 
returned to mass storage. 

GRAPH subroutines may be divided into three 
categoriss: 

1. DRIVER--provides for initialization of the progran 
and drives subroutines. It also passes control to 


the next program in the OS. 
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2. AXIS and Labling Subroutines--establishes the axis of 
all graphs and lables as appropriate. Provides for 
formatting of descriptive parameters, aS appropriate 
builds grids for graphics. Six subroutines comprise 
ES.  GrEOup. mney are SALTS, TICKS, AXLAB, TITLE, 

GRPAR, GRID. 

3. Inverse Scattering--Subroutines in this section graph 
the results of processing in INPUT and MATH. They 
reproduce the four basic input waveforms, the Fourier 
transforms of the Direct, Backscatter and Impulse 
Response, the Physical Optics Shape of the target, 
and the time domain ramp response. Four subroutines 
comprise this group. They are GRIN, GRFOUR, GROPP, 
and GRAMP. 

A general flow diagram for GRAPH is provided in 
megure 4.9. A disscussion of all Subroutines in GRAPH 
follows. 

a. Graph DRIVER 

As for INPUT and MATH DRIVER, GRAPH DRIVER is 
meme Draneipel initializer and contzoling routine for this 
program, coordinating the flow of processing. Unlike MATH 
DRIVER, GRAPH DRIVER does not of itself perform any data 
processing. That function is retained for the subroutines 


exclusively. 
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GRAPH DRIVER is Located in lines 100-1190 and 
Comprises 19% of the total progran. GRAPH is completely 
independant of data input from the other programs in the OS. 
All inputs are from mass storage Tape 2 or the GS keyboard. 
All results are returned to CRT display of the 4052. This 
routine is highly operator-procsssor interactive. Appendix 
F describes a typical run with messages as produced. 

GRAPH DRIVER is always initiated in line 100 
when being entered external to the GRAPH Program. It is 
initiated from line 290 by returns from internal 
subroutines. Lines 100-300 initialize the program and input 
parameter data from mass storage. Lines 320-750 develope 
the correct Menu, With the function/subroutine selected 
input to PO in line 760. 

The Henu that is used in GRAPH Tis shown in 
Apendix f°. Most of the fiunctions on the menu are 
self-explanatory. A few remarks concerning items 10, 11, 
12, and 13 should suffice to rémove any ambiguity. 

hecge 10, “CONTINUE OrTigtza lion", initiates the 
Riad low-pass filtering optimization of the impulse response 
in lines 950-1030. Program control with the operator 


designated optimization factor is input to Program MATH. 
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Item 17, "NEXT TARGET", passes program control 


Semeehe ACjUgscation of ah entirely new target. 


Fh 


so. LNPUT 
Parameters are retained and passed along with prograa 
Gentrol. 

TEem 12, TEEN LSe oD GRAPHING", terminates 
processing. The program can be reinitiated by typing "RON". 

Lines 820~900 select t*+he appropriate subroutine 
and funetions. Line 940 returns to the menu for further 
selection for all items except 10,11, 12, 13. 

The next series of subroutines to be discussed 
Betcroragithe graphics axas formatting, labling and titling. 
These subroutines set the environment tha* makes the 
information of use td an observer. POLS PEGKS ,- AXLAS, wand 
TITLES are subroutines based on programs vreviously developed 
ye TeXtronixy, Inc. fRef. 22j] All have been modified toa 
greater or lesser degree dependant on the needs of the TDRL. 

Dee SubEOUtiInes AXIS/TICKS/AXLAB 

SUp=oucine AKIS drives Subroutines TICKS and 
OAS . They perfora the general axis definition function 
BOoeeene Output graphics. Subroutine AXIS serves as the 
imput portal for data from the function subroutine that was 
chosen by the operator. The results returned to the CRT of 


the GS are "neat" ticks, labled and scaled to the proper 





Values of the input function. Subroutine AXIS occupies 
Statements 2000-2250, TICKS is in lines 2300-2610 and AXLAB 
2S in statements 2700-3190. These three subroutines occupy 
18% of program space. 

Data input from the function to be grapked are 
the minimuna and maximum horizontal values of the viewport in 
Yt and V2, and the mininum and maximum vertical values of 
the viewport in V3 and V4. These constants are in Graphic 
Display Units (GDU). GDU's are a measure of the resolution 
capability of the SS screen. The four arguments defined 
above refer to the actual limits for the graphic information 
Seetne display. The GS will draw no lines outside these 
eamits. 

Also input to Subrowtine AXIS frome the driving 
function subroutine are the ainimum and maximum horizontal 
values of the function itself, wW1 and #2, and the ninimun 
and maximum vertical values, W3 and W4&. These values are in 
User Data Units (UDU"'s). hiat 1S; in units that the user 
defines, such as nillivolts, seconds, meters, etc. 

Pinally, AXIS receives the value of the flag SO. 
This specifies whether the vertical axis (SO0=0) or 
horizontal axis (S0=1) is being labled. The magnitude of the 


Menrezontas Scaling Zunction for the input wave, S1 (received 





from INPUT), and the definitions 


axis and horizontal axis are required. 


The performance of AXIS, 


adequately described in the 


(Ref. 23] This repoct will describe cnly 


which have been changed or modified to 
Sucput for TDRL purposes. 
Principally, the 


Which the axis is lable with UDU's 


These changes lie in statements 2860-3190. 
must be checked: 


1. Horizontal (S0=0) 


TEKS, 


Tektronix sources 


changes effect the 


see Uspereiheakera 


Or vertical axis (S0=1) 


in JDU's £or the vertical 


and aAXLAB is 
referances. 


those statements 


produce the desired 


manner in 


AXLAB. 


Three conditions 


lables? 


2. Graphics of a Fourier transform? (4<P0<8) 


3. Check amplitude of maxinun 


UDU's. 
TABLE IV defines the re) 
Statements in Subroutine AXLAB and 


grapned. 
Cc. Subroutine TITLE 


This subroutine positions 


graphic above the plotted values. bee 1S 


3200-3310 and comprises 2% of the progran. 
and the 


maximum norizontal coordinate wW2 
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vertical value 


the 


(P6) in 


ationship between 


the functions t0 be 


title of the 
located in lines 


The value of the 


graphic computed 
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Graphic Labling Relationships--GRAPH AXLAB 





























| GRAPH mi g aay REMARKS _ 
| { 
GRIN HORIZONTAL ; HORIZONTAL | 
| GRAM AXTS AXIS : 
{ 28 80-2890 $1*1021000. Lable Axis 
| 29 10 $1*10<1000. Ladle Axis 
| | 3090 { Position Dimension Lable | 
| | 3110 Dimension Lable | 
roam a SSS ee a ee ey 
3140-3160 Print Dimension Lable a 
i i 3180-3190 | Return | 
| --+-—_ -—_------ OO 
| 
| 
| 
i 
i 
















Li = = 
VERTICAL VERTICAL 
AXIS i AXIS | 
- Start | 
2940-2970 P(6)<1. Lable axis | 
-~-—-—-—_-—__-—___- +- —_-- SS Se 
| 2980-2990 | 1$P(6)<10. Lable axis | 
| 3000-3010 i: P({6)210. Lable axis | 
} 3130-3150 | Print Dimension Lable 
2 eee ee ee 
| 3180-3190 | Return | 
HORIZONTAL «HORIZONTAL 
AXIS AXIS 
| 2861-28562 t 
| 2863-2864 "Magnitude" Axis Label = 
| -——— | -—- ———- ——-- -—-- eee 
| {| 2865-2866 | hase” Axis Lable 
| [— 3171 i Print Dimension Lable — o 
| | 3180-3190 } =a i .~~C~CS«~<a 
Voge LC Ls VERTICLE © 
{ AXIS \ AXIS 


Return 


a 8 OOO eet 
| 2860-3170 —— Same as for Horizontal Axis | 
Pf 3173-3190 Pe Retarn 


NOTE: All graphic routines use lines 3020-3040 
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vertical maximum P(6) are input by the graphic subroutine. 


ie | 


The graph title centered over the graph is returned to the 


CRT. 
d. Subroutine GRPAR 
GRPAR graphs the parameter input by INPUT and 
developed by MATH. This subroutine is located in lines 


3400-3870 and comprises 8% of the progran. All graphs have 
the run number, date, target distance from transmitter, 
antenna number, target number, averages taken, optimization 
and any specific remarks printed. The specialized 
parameters printed for particular graphs alone with the 


statement number apolicable are given in TABLE VY. 


TABLE V 


Specialized Parameter Prints--GRAPH GRPAR 


GRAPHIC | - PARAMETERS | LINES ~ 





SED EE ES DRTC EDAD atte ent oR PU SNEED at ee! 
| Gary [MAXIMUM PEAK VALUE sditaa 
| caamp [ HINTMOM PEAK VALUE | | 
| (ne | 
| MEAN VALUE ; | 
————————— 


7 Gaee t DIAMETER (INCHES) Rs) eee So 


el 
| ia WIDTH (INCHES) © | | 


150 








2. Subroutine GRID 

This simple subroutine draws referencing grid 
lines on the appropriate graphic. The computed graph 
maximum and minimum are input from the graphic subroutine 
mong With the axis tick intervals. Grid lines are output 
eo the CRY from GRID. Lines 3900-4010 Contain wene 
Subroutine. GRID requires 2% of program space. 

[hs concludes the iscussion of general 
Meapoics formating subroutines. Ths @r2ving graphic 
Subroutines will now be discussed. 

Subroutines GRIN, GRFOUR, and GRAMP are all 
Similar in formatting structure. With some judicious choice 
of progranming flags and routine definition, it should not 
be too difficult ts combine these three subroutines into a 
Single driving graphic subroutine. GROP is formatted 
differently, but still retains a few features common to all 
the graphic driving subroutines. The similarities between 
the subroutines will be discussed first. The carticular 
differences applicable will be discussed secondly. 

Each graphic subroutine is composed of four 
basic blocks. They are: 

1. Input waveform data from mass btorage 


2. Define waveform and graphic boundaries 
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3. Draw the graph of the function using special defining 
functions when necessary 
iomectolane Che Graph with the correct title(s) 

GRAPH DRIVER inputs the function value PO which 
is used to identify the location of the data storage file in 
which the appropriate data rests. This single row array is 
usually returned in X. The waveform window and viewport 
parameters, bounding the output graphic display, are next 
determined. The maximum horizontal axis values are found by 
Mutiplying the distance between a single data point (in 
UDU's) by the total number of data points in the array along 
the X-axis. For those waveforms listing the mean graph 
values as part of the descriptive parameters, the sum of all 
the array element values is obtained and then divided by the 


discrete number of array points: 


tmam 
— 


x(1)/N Cai 


To ensure that the function is graphed along the zero axis 
of the formatted graphic, the means are subtracted from the 
array elements. In similar fashion, the root mean square is 


found: 
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N 
RMS = W6 = [ £ x(i)*/512]1/2 (4.12) 


1=1 


The maxinum and minimum values for the input 
data are found by implementing the 4052 call "MAX" and "MIN* 
Boucines. These values are returned in WwW4 and wW3 
respectively. The values are then output to the axis 
defining subroutines for e2stablishment of appropriately 
scaled and labled axis. 

Next, the GS cursor is positioned to the 
beginning of the array with a "MOVE" command. For the 
humber of discrete data in the array, the function is 
graphed on the GS Cart. Pinally, the graphic is titled and 
identified by the appropriate parameters. Figure 4.10 
PaWihdes a general flow diagram of the graphic driver 
routines. 

The particular aspects of each graph will now be 
briefly discussed. 

ff. Subroutine Grin 

This subroutine is the driving graphic for the 
display of the inputted Direct, Incident and Augmented 
Waveforms, and derived Backscatter waveforn. It resides la 


Statements 5000-5410 and makes-up 7% of the progran. 
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Dimes general descriptiom of graphic routines 
provided above is most typified by GRIN. 

TwO special notes need to be made. berst, 
graphing begins with the ninth element of the array. This 
is to prevent the display of an erroneous spike input by the 
DPO. (These erroneous values have a minimal effect on the 
processing). Secondly, the Function is graphed in discrete 


steps equal to 1/N of the total array value along the axis: 


Graphed Point = (Q * 21), X(Q) C4, 13) 


where: Q array number 


X (Q) array element 

The title of the graph is selected by the 
routine located in lines 5290-5380. 

The horizontal dimenion for the graphic is time 
while the vertical axis is scaled in volts. 

ge Subroutine GRFOUR 

This subroutine is also fairly typical of the 
general graphic method. It displays the magnitude and phase 
components of the spectral data in polar forn for the Direct 


and Backscattered waveforms and for the impulse response of 


the target. 





Only the first 60 elements of the nagnitude 
array are graphed. Values beyond the first 60 are ofa 
random nature due to the processing performed to derive the 
transforms. Spectral content much beyond 10 GHZ is, 
therefore, negligible. All values of the phase array are 
plotted, although the note concerning the randomness of the 
data is also applicable here. 

GRFOUR is actually performed in two distinct 
operations. The first graphs the magnitude of the 
transfora. The horizontal axis dimensions are converted to 
frequency from time by Finding the frequency of each data 


BoOlnt : 
Graphed Point = Q/(Z1¥*257), X(Q) C4. say) 


The second operation graphs the transform phase. 


The horizontal axis dimensions are converted to radians: 
Graphed Point = (Q*27¥)/257, X(Q) CAS) 


GRFOUR is programmed in lines 5500-6080. It is 10% of the 
eOtal program. 
h. GROP 
GROP is the graphic that reproduces the shape 


representing the target. That is, the target shape based on 
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Physical optics principals and uncorrected for field effects 
induced by residual currents. ite represents a 
two-dimensional cross-sectional area cut of the target. 
Therefore, only target length and height are shown. 

This subroutine is constructed generally in the 
same format as the other graphic subroutines. One important 
difference is that a true one-to-one relationship must be 
established for the coordinate axis. However, before this 
is done, the netric dimensions of the target length must be 
bounded. This is done in lines 6210-6240. Note that the 
array BO, source for the shapes storage in MATH, has already 
been scaled for the metric dimension in the y-direction or 
amplitude. The x-dimension, or length, is scaled for metric 
measures in line 6220. This is accomplished by multiplying 
the time between samples by the total number of samples in 
the array representing the target. This value is then 
compensated by the averaging antenna parameter established 


in MATH, SS: 


10 
Merce) 1.2 
we = CODE) rr exto (SS) (4.16) 


where: 2Z = sampling time between single samples 


Sah 
tl 


measured Length in number of samples 


S5 antenna length compensator 
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Q = low-pass filter compensator 

Q is always one for zero optimization, or no low-pass 
me tering. It has effect only when optimization is 
performed, and is designed to compensate for the spreading 
of the target dinensions in the X-darectron durang 
filtering. Note also, that all defined targets are measured 
in inches. All outputs are scaled inches. + should be a 
relatively simple procedure to generalize the program to 
handle either metric or english units. 

The next najor difference between GROP and other 
graphic routines is that all graphing is done in GROP itself 
Without calling AXIS or any other external formatting 
Subroutine. As mentioned above, the important concept is to 
properly relate the co-ordinant axis ona ‘true one-to-one 
Blot. This is done is statements 6290-6360 and is called 
"axis scaling”. Axis scaling is done in two steps. FLOSt 
the general x-axis scaling factor AO is found. This is done 
in 6290 by dividing the maxinum horizontal bounding UDU by 


the viewport established by the display: 


AQ = (W2 - W1) / (V2-V1-5) (41.41) 
This yields UDU's/GDO inthe x-direction. The y-axis 
Sear ng cactor Al is then fixed to this value. The 


tBexs) 








dimensions in the y-direction are then checked by the 
program in lines 6320 to ensure that "clipping" or windowing 
out of the data does not occur. This is necessary as the 
vertical viewing axis is only 62% that of the x-axis. If the 
y-dimension were more than 52% of the x-dimension, data 
would be lost during graphing. To prevent this, the x-axis 
values are fixed so that the maximum x-dimension will always 
be plotted. This is an effective reduction in the size of 
the total displayed shape. To retain the ‘true one-to-one 
axis relationship, the y-direction dimensions must be 
reduced by a corresponding percentage. This reduction 
factor is computed in line 6330 and applied to the array in 
644d, Note that in 6330 twice the maximum y-dimension is 
used. This is done because the graphic output produces the 
Mirror image of a axisymmetric object while the array itself 
1s single-sided. 

Once scaling is established in statement 6380, 
the array is plotted. As data is graphed it must be 
dimensioned in english units, compensated for optimization 
effects, and compensated for any viewport clipping. This is 
done on a point-by-point basis in line 6410. Note that the 
acray has already been totally dimensioned and compensated. 


Tt is in fact the plotting point which now must be taken 


| Sie, 





into consideration. This compensated plotting point is 
returned in 01 and is used in drawing the graph. ginaily, 
the mirror image of each point is graphed to produce a 
"true" image. A flow diagram for GROP is given in Figure 
ott. 
i. Subroutine GRAMP 

GRAMP graphs the time domain ramp response in 
lines 6600-6910. It comprises 54 of the program. 

GRAMP is identical to GRIN in practically every 
detail. Note that the array X is multiplied by 1000 in line 
6670. This converts array amplitudes to millivolts. 


GRIN also employs this conversion. 
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ieee nis chapter, the results of several selected 
measurements will be presented. The observations made were 
chosen to demonstrate the flexibility of the Operating 
System as well as verify experimental results with the 


theoretical. 


A. MEASUREMENT 1--VARY BORESIGHTED TARGET DISTANCE 
1. Deseription of Setup 
The purpose of this measurement is to observe what 
effect locating the target a distance different from the 
transmitting antenna then that specified for the receiving 
antenna, but maintaning the boresighted relationship of the 
antennas has on measurements taken by the OS. The 
receiving antenna is placed 2.54M from the transmitting 
antenna. The target is defined from the target library on 
Tape 2. 
2. Results 
The results are given in Figures 5.1 to 5.3 with the 


backscattered waveforn, ramp response and physical optics 


Shapes being presented. 








The ramp response is seen to be typical for a 
cylinder although somewhat reduced in size. This is to be 
expected since the propagation attenuation losses are 
greater for the distance at which the target is located in 
this experiment than it is for the normal control case. 
With the distance taken into consideration by the Operating 
System software, it can be seen that the resultant physical 
optics shape is a close approximation. The target width 
(length) is within 10% of the true length and the diameter 
is 15.4% of the true diameter. The uncorrected width is 49% 
greater than the true length. The backscattered waveforn 
Shows the effect of some spreading in comparison to the 
control case and also displays a reduced S/N. 

3. Conclusions 

The distance between the receiving antenna and 
target can be arbitrary andthe Operating System will be 
able to provide a reasonably accurate interpretation of the 
Baiada . As the distance increases, this accuracy can be 
axpected t9 decrease directly to the limit of the Systems 


ability to seperate the target and noise. 


Tee, 





TARGET: RUN DATE DIST ANT TGT REMARKS 
2 9-7-8! 2.544% 1 i LOCATION DIFFERENCE 


MAXIMUM PEAK UALS coc asc sccece 41,68 wvV 


MINIMUM PEAK VALUE... cece ccees “1.49 MU 
RMS VALUE. cer eccccscssccssvee +8.62 mV 
MEAN VALUE... cccccsvccccvevces @.98 a 
NUMBER OF WAVEFORMS AVERAGED = 52 OPTIMIZATION VALUE # 8 





8.58 1.00 1.509 2.00 2.58 3.88 3.50 4,88 4.59 3.88 nS 


Pigure 5.1. MEASUREMENT 1--Backscattered Waveform. 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
RSE) 2.54m 1 1 LOCATION DIFFERENCE 


MAXIMUM PEAK Ue eee +6.86 mU- 
MINIMUM PEAK WEGIBIE cietae esies «66 0 -196.33 mU-N 
RMS UiGiimliite cis'etete coo ees eeeeeeees 112.76 a Od | 


MEAN Ua soc bee ees 6 See eee -90.84 mV-t4 
NUMBER OF WAVEFORMS AVERAGED = 52 OPTIMIZATION VALUE = @ 





8.58 1.08 1.532 2.08 2.58 3.08 3.58 4.800 4.58 5.38 nS 


Pigure 5.2. MEASUREMENT 1--Ramp Response 
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TARGET: RUN eae DIST ANT TGT REMARKS 
eae aot Dy ee | t LOCATION DIFFERENCE 


OTAMETER <INCHES=-CORRECTED) 6.9525 ¢ 6.89890) 
WIDTH CINCHES~CORRECTED) 13.2847 (12,8880) 
WIDTH CINCHES~UNCORRECTED)&x% 17.3711 (12,0880) 
NUMBER OF HAVEFORMS AVERAGED # ay OPTIMIZATION VALUE =» 9 


RAMP RESPONSE <= “PHYSICAL OPTICS 


Pigure 5.3. MEASUREMENT 1--Physical Optics Shape 
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Bo. MEASUREMENT 2--TARGET AND RECEIVING ANTENNA NOT 
BORESIGHTED 


1. Descriptio 


of Setup 

An attempt 1s made to observe the effect of placing 
a target at an arbitrary location on the [DRL ground plane. 
The transmitted field is not axially incident on the target. 

For this experiment, the target (cylinder 1) was 
placed 2.11 meters from the transmitting antenna and ata 
909 angle (approximately) to the receiving/transmitting 
antenna axis. The receiving antenna is located in the 
control area of the TDRL. The target is at about a 600 
angle to the receiving antenna and located 2.44 meters from 
ant 6 The experimental setup is diagramed in Figure 5.4. 

2. Results 

The results ar2 given in PaeGguiies 529 CO S67 
representing the ramp response, physical optics shape and 
backscattered waveform. Note that the distance input to the 
systam is incorrect. The correct value as read from the DPO 
should be 1.57 meters. The diameter should be multiplied by 
meme factor 1.53. Doing so yields a diameter of 6.3341 
inches as the true corrected value. The resulting 
compensated values are thus very close to the targets actual 


dimensions. These figures may be suspect as subsequant runs 
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TRANSMITTING 


1-2) Mm 


RECEIVING 
ANTENNA 





Figure 5.4. Measurement 2--Physical Setup 


never yielded the same degree of accuracy, although they 
approached it. The important item is that the physical 
Optics shape is a close approximation of that expected to be 
representative of a cylinder. 


3. Con¢elusions 





It would appear that the targets location can be 
generalized if the shape only and not the dimensions is the 
seonitcacant item of interest. However, time precluded 
further verification of this observation. It cannot be said 
with absolute certainty that this observation is completely 


general. Other configurations with the target located at 
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different distances from the boresighted transmitting and 
receiving antennas should be observed. Note that for proper 
compensation for propagation attenuation, the distance the 
target is from the receiving antenna is determined by 
reading the DPO Distance Scale. This should be the 0S 
distance input value vice the straight line measured 


distance from transmitting antenna to target. 
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TARGET: a“ agTe DIST ANT a REMARKS 


7-81 1.27% 1 NOT BORESIGHTED 
MAXIMUM PEAK VALUE... cc ccc wane +47.50 mU-M 
MINIMUM PEAK VALUE. ec cece scans -329.05 mV-¥ 
RAS VALUE. .cccccccccccscccvees 181.31 mU-M 
MEAN VALUE. ccccrccccsccsccccce “133.00 mU- 
NUMBER OF WAVEFORMS AVERAGED = 4 OPTIMIZATION VALUE = @ 






cal ho 


8.58 1.98 1.539 2.8008 2.58 3.908 3.58 4.00 4.50 3.98 ns 


Pigure 5.5. Measurement 2--Ramp Response 
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TARGET: RUN DATE Ost ANT TGT REMARKS Tite ne 
3 98-7-81 teers 1 L NOT BORESIGHTED . 
DIAMETER CINCHES-CORRECTED) 4.1447 ¢ 6.9898? ae 
WIDTH CINCHES-CORRECTED> 11.8218 (12.9880) ; 
WIDTH <INCHES-UNCORRECTED) &% 15.8789 (12.3000? 
UMNBER OF WAVEFORMS AVERAGED = 4 OPTIMIZATION VALUE #= 86 





Pigure 5.6. Measurement 2--Physical Optics Shape 
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TARGET: RUN DATE OIST ANT TGT REMNARKS 
alt dead - P| sor 4 1 


1.2 NOT BORESIGHTED 
MAXIMUM PEAK VALUE. cccccccnsee 42.48 mV 
MINIMUM PEAK VALUE... cc accccves 2.439 mY 
RMS VALUE. .cccccccccccvcsccene 48.98 mV 
MEAN VALUE... cece eone eeccceeen 8.88 mV 
NUMBER OF WAVEFORMS AVERAGED = 4 OPTIMIZATION VALUE = @ 


BACK SCATTER WAVE FORM 





2 | a a a a a a 


9.58 1.08 1.58 2.00 2.58 3.80 3.58 4.88 4.58 5.98 ns 


Figure 5.7. Measurement 2--Backscattered Waveform 
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Ge MEASUREMENT 3--ANTENNA EFFECTS ON TARGET RESPONSES 
1. Description of Setup 

Thre purpose of this neasurement is to verify that 
the effects of the antenna are in fact deconvolved out 
during processing. “Most previous runs have been performed 
uSing antenna 1 (TEM Horn). This measurenent is performed 
using three different length monopole antennas. The 
antennas used were antenna 2 (see TABLE ITI) (85 mn 
monopole), antenna 8 (205 mm monopole) and antenna 12 (285 
mm monopole). The target used was target 1 (see TABLE I). 
Both target and antenna were placed in the control 
locations, 1.27 meters from the transmitting antenna. 
21-acquisitions were averaged for each run. The results 
obtained should be compared with those obtained previously 
mor the TEM horn. 

2. Results 

The results are presented in Figures 5.7 to 5.24 
Each antenna has the results of its run shown in six graphs. 
The first graph depicts the direct waveform as seen by the 
individual antenna, the second 1s the ramp response 
resulting, the thivecd) f2guce 2s the physical optics 
representation of the target shape, and the final three 


figures are the optimized fraquency domain impulse, time 
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domain impulse and ramp response for the target. Figure 5.7 
BO 5.12 are for monopole antenna 2. Puguresmeo to tous. to 
are for monopole antenna 8 and figures 5.19 to 5.24 are for 
monopole antenna 12. 

The first obvious result is that in comparison of 
the ramp responses, the smaller the antenna is in relation 
to the target, the less typical is the ramp and, 
correspondingly, the physical optics shape. However, for 
antennas considerably taller than the height of the target, 
a Significant amount of noise pollutes the ramp response. 
In comparing the target impulse responses, it can be seen 
that the exact opposite appears to be true. The most 
representative impulse is for the the shortest monopole. 
The frequency domain impulse graphs show that the longer 
monopoles accept a higer frequency content than the shortest 
monopole. 

3. Conelusions 

Euas t, he closer the antenna is to natural 
resonance with the incident pulse, the more faithful appears 
the resulting impulse response to theory. But it does not 
naturally follow that the corresponding ramp response will 
be as true a repsentation. The opposite appears to be se 


It may be that the short antenna is shadowed by a portion of 
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the target. More likely, the short antenna is just not 


large enough to accept the full radiated backscatter of the 


target. Conversely, the 2385-mm monopole is sensitive to 
environmental noise. As can be seen in figures 5.20 and 
oe 14, the ramp responses display the relatively greater 


Sensitivity of the 285-mm monopole to environmental noise 
than the 205-mm monopole. Although the observation is 
limited to just these runs, it may be possibly concluded 
that there is a "best" monopole antenna length for a 
particular target size. 

In conclusion, it can be said that the effects of 
the antenna are in general deconvolved out. However, other 
effects iid rectly related to the target/antenna 
configuration may introduce factors which materially effect 
the reproduction of the target shape and its overall 


responses. 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Weweaii=81 t.2°7M 2 1 DIPOLE ANTENNA #2¢85MM) 


MANIMUM PERK VALUE... coeccsces +76.01 mV 
MINIMUM PEAK VALUE... .secccees “91.33 mV 
BOUL UEs soc c sce o cess ccvcsees +#31.83 mY 
MEAN VALUE. ccc ccccccscvcccccces 8.00 mV 
NUMBER OF HAVEFORMNS AVERAGED 2 21 OPTIMIZATION VALUE = @ 


GIRECT WAVEFORM 


ROLLE 
TCA 
OMT 
PAV CCE eee 


9.50 1.08 1.50 2.88 2.58 3.90 3.58 4.88 4.58 5.938 nS 
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Pigure 5.8. Measurement 3--Direct Waveform (85 mm) 
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TARGETS RUN DATE OIST ANT TGT REMARKS 
eer -8l 1.27n 2 1 DIPOLE ANTENNA #2¢85MMD 


MAXIMUM PERK VALUE... ccc eeeeee 782.43 MUM 
MINIMUM PERK VALUE... . ec cceeee 7909.93 MUM 
RMS VALUE. cc ccc c er cccccceccces 3742.88 mUnM 
MERN URLUE.. .ccccccescccececes *153.01 mU-M 


NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 6 


RAMP RESPONSE -- TIME 








Re 
oo] | rey | 


8.58 1.08 1.50 2.808 2.58 3.88 3.58 4.88 4.58 3.88 ns 





Figure 5.9. Measurement 3--Ramp Response (85 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
fo catte=ot leer 2 | DIPOLE ANTENNA #2¢85MNMD 


DIAMETER CINCHES-CORRECTED) 6.6879 ¢ 6.8989) 
WIDTH CINCHES-CORRECTED) 12.5283 (12.8008) 
WIDTH CINCHES-UNCORRECTED) «x 10.1367 (12.0608) 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 


RAMP RESPONSE -- PHYSICAL OPTICS 


f ( 


Sha 
| eet 


Pigure 5.10. Measurement 3--Physical Optics Shape (85 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
WeooUl=Sloel.2rn fe 1 DIPOLE ANTENNA #2¢85MND 


MAXIMUM PEAK VALUE.....cccecee 9.90898 mV 
MINIMUM PEAK WRU crete ccc 6 ees @ -8.0007 mY 
RMS JRL Ar 8.8982 aU 


MEAN VALUE. «- 8.9000 mV 
HUMBER OF WAVEFORMS. AVERAGED 2 21 OPTINIZATION VALUE = 8.5 


2 OR ap o> GP > 4 Ge OD 6 SP 4 o® 6m @® 6m 6D O® OD S® 4 a> OD 4 4b O® 4 4D 4 OD OD 4 4D OD GE OD OD 4 OB OD om 4 4 4 4 o> 6B 6D OD 6D OD OD OD om OD OD OD 42 SE OD &® OD o® OD OD GF O@ GD EB ED 


8.8885 





8.8888 






PUTTIN 
SSR 


@.58 1.08 1.58 2.08 2.58 3.38 3.58 4.08 4.59 5.98 nS 


-8. 8885 


surenent 3-~-Tine-Domain Iapulse Response 


Figure 5.11. Mea 
~ (85 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Weme-1i-8i 1.2°h 2 | OIPOLE ANTENNA #2¢85MM)> 
NUMBER OF WAUEFORNS AVERAGED 2 21 OPTIMIZATION VALUE # @.5 


SSS BSW SS SOs GSO SOP Oe 22 2 22 2822 2228 22 O22 2228 PSV 2 O82 2B 2 OOO DWlA ®OOBS_S*Q_ LD B2@ @ @ @ @ @ @& o& oF 4® a =o 


IMPULSE FFT--MAGNITUDE 


0.840 
9.920 


8.869 


8,08 2.00 4,09 6.98 8.08 18.08 12.88 GH2 


TELL 
Cl cell 
+P ELLE etter 


8.8 8.5 1.9 1.5 2.8 2.5 3.9 3.5 4.8 4.5 5.9 3.5 6.8 RAD 









8.88 





-5. go 


Pigure 5.12. Measurement 3--Prequency-Domaain Impulse Res- 
ponse (85 mm) 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
Poesatl=Si 1.27% e 1 DIPOLE ANTENNA #2¢85MMD 


MAXIMUM PEAK VALUE... csccceceee *613.71 MU-M 
MINIMUM PERK VALUE... ce cece eee = 7395.20 AUR 
RMS VALUE. cc coer cc cresevecccee 284.65 MU-H 
MEAN VALUE... ccc crcccccsccsesee 141.780 MUM 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 2.5 


RAMP RESPONSE -- TIME 


9 a a 
TCC CCAZINE 
8 
PEC CEE CCRT 
per et 
Tw itt ttt lo 


0.59 1.08 1.58 2.98 2.58 3.88 3.50 4.90 4.58 3.88 nS 
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Figure 5.13. Measurement 3--Ramp Optimized (85 mam) 
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TARGET: RUN DATE OfST ANT TGT REMARKS 
feee=ii=st i.27h 8 t DIPOLE ANTENNA #8 (205MM) 


MAXIMUM PERK VALUE... ..cccevee +104.92 MU 
MINIMUM PEAK VALUE... cee ceeeee TIL6.61 mV 
RMS VALUE. coc ccrccvavcccvcscscce +59.58 mV 
eH VUE. oc cc etc c cs oascce> 3.80 mU 
NUHBER OF WAVEFORMS AVERAGED 2 21 OPTIMIZATION VALUE = @ 


DIRECT WAVEFORM 


eee | 
JT 
eo 7) TN 
| fF | | | 
Nene 


9.59 1.88 1.58 2.80 2.58 3.00 3.50 4.808 4.58 5.80 nS 


aV 








igure 5.14. Measurement 3--Direct Waveform (205 am) 
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TARGET: RUN DATE OST ANT TGT REMARK 


Memeo i=-Siva.2rh 8 
MRAIMUM PEAK VALUE... cccccee 
MINIMUM PEAK VALUE... ccccsccves 
RTS PRUE. 6 wc occ cece asc vscvee 
Mery VRE. . csc cc ects ec ccene 


NUMBER OF WAVEFORMS AVERAGED = 21 


1 DIPOLE ANTENNA #8 (285MM) 
+5.83 mU-M 
-609.89 mU-N 
#357.09 mU-M 
-299.54 mU-M 
OPTIMIZATION VALUE 2 @ 





9.58 1.8908 1.58 2.00 2.58 3.88 3.58 4.30 4.58 3.98 nS 


mpm ee me ee 


Pigure 5.15. Mea 


surement 3-—~Ramp Response (20 


5 mn) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
V4e S-ii=Si i.27M 8 t DIPOLE NEA @8 (295MM) 


OITAMETER CINCHES-CORRECTED) 3.3341 ¢ 6.0000 
WIDTH CINCHES-CORRECTED> 12.0734 (12.9098) 
WIDTH CINCHES-UNCORRECTED) 4% 19.3945 (12.6008) 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 


RAMP RESPONSE ~~ PHYSICAL OPTICS 








Figure 5.16. Measurement 3--Physical Optics Shape (205 an) 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
Meacoii=81 (.27M 8 ! DIPOLE ANTENNA #8 (205MM) 


MAXIMUM PEAK VALUE... ccccceess 9.90806 mV 
MINIMUM PEAK VALUE.......-5006 “8.0806 mV 
POTERURIMUEGi sc csesseeccccsces+s 0.9002 mV 
MEAN VALUE..... 8.9000 AV 


NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8.5 





8.58 1.830 1.50 2.80 2.59 3.00 3.58 4.86 4.58 5.88 ns 


Figure 5.17. ' 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Demetri i-3) 1.27% & { DIPOLE ANTENNA #8 (265MM) 
NUMBER OF WAYEFORMS AVERAGED *® 23 OPTIMIZATION VALUE ® 6.5 





8.09 2.898 4.00 6.98 8.08 10.98 12.88 GH2 





8.9 9.5 1.9 1.5 2.8 2.5 3.8 3.5 4.9 4.5 5.9 5.5 6.8 RAD 


Figure 5.18. Measurement 3--Prequency-Domain Impulse Res- 
ponse (205 mn) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
14 8-11-91 1.276 8 L OIPOLE ANTENNA 8&8 (285MM) 


MAXIMUM PEAK UALUE.......0062- #0419 MU-M 
MINIMUM PEAK UALUE........-++. -537.10 mU-M 
RMS WALUE.....cescceccvecccves +340.09 nU-M 
MEANMUGEUE..<..0cecseocesceses -25503! nU-H 

NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 2.5 





@.59 1.00 1.50 2.88 2.580 3.80 3.359 4.90 4.50 3.388 nS 


Pigure 5.19. Measurement 3--Ramp Optimized (205 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Resor i=St i.27M i2 1 DIPOLE ANTENNA #12¢285MM) 


MANIMUM PERK VALUE. cc coc ecceee§ *101.31 mV 
MINIMUM PEAK VALUE... cece ecceee 7122.78 mV 
RMS JALUE. . ccc cc ccccevesenves 67.52 mV 
MER VALUE. . cccccccvccsccceses 9.80 mV 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = @ 


OIRECT WAVEFORM 


TT TTT 
TSE 
“J a 
2 
COTE 
es 


0.50 1.88 1.580 2.08 2.58 3.090 3.30 4.90 4.58 3.80 nS 
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Figure 5.20. Measurement 3--Direct Waveform (285 mm) 


188 





TARGET: RUN DATE OIST ANT TGT REMARKS 
fems-il=31) l.crm te | DIPOLE ANTENNA #12¢285hM) 


MAXINUM PEAK VALUE... cc cccecees #£39,.35 mVU-M 
MINIMUM PEAK VALUE... ccc cece ee “SOSe13 MUM 
RMS URLUE. core cc cccccrccvcenee 320.92 NUM 
MEAN VALUE. .cccccccccccrccsees Sldcerea MUaH 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = B 





@.56 1.98 1.509 2.08 2.58 3.88 3.58 4.99 4.50 5.89 nS 


Pigure 5.21. Measurement j3--Ramp Response (285 ma) 
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TARGET: RUN SATE OIST ANT TGT REMARKS 
feee-li=31 t.27M 12 i DIPOLE ANTENNA #12¢28S5MM) 


DIAMETER CIHNCHES-CORRECTED) 6.0000 ¢ 6.88980) 
WIDTH CINCHES-CORRECTED> 12.8000 (12.9888) 
WIOTH CINCHES-UNCORRECTED) *% 16.2395 (12.9080) 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 


RAMP RESPONSE -- PHYSICAL OPTICS 
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Pigure 5.22. Measurement 3--Physical Optics Shape (285 mom) 
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(y: 


TARGET: RUN BATE OIST ANT TGT REMARKS 
Pemeairi=siisert t2 1 DIPOLE ANTENNA 812¢285MM) 


MAXIMUM PEAK VALUE... ccc cece ee 8. 00806 MU 
MINIMUM PERK MSIE ooo ei oo oleceue evece -8,0803 mV 
RMS MAbs) ereletetctergyesicte'e 4s 6 66 ele 6 60 9.38882 mV 


ERM PORLUE ce cco ene cseeeteece 8.08090 mU 
NUMBER OF WAVEFORMS AVERAGED 2 2 OPTIMIZATION VALUE = 6.5 


es 6 


IMPULSE RESPONSE 


AERO 
ACCEL 
» LSMARABeM 
JERSE 


9.58 1.00 1.50 2.08 2.589 3.08 3.58 4.80 4.58 5.880 nS 





8.0985 





8.8898 









-8.98835 





asurement 3--Time-Domain Impulse Response 


} ® 3. Me 
Figure 5.2 (335 mm) 
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TARGET: RUN ae O1ST ANT TGT REMARKS 
femoetsoet) 6 6ll. 27m 6 612 ft DIPOLE ANTENNA #12¢285HM) 
NUMBER OF WAVEFORMS AVERAGED 2 21 OPTINIZATION VALUE = @.5 





5.80 2.08 4.08 6.26 3.38 10.08 12.868 GHZ 


INPUCSE ERT==RHAce 


eR REE 
Pet acl hat Uh 
TT YT pana OT TT 
tee 


G50, 6.5 1.0 1.5 2.9 2.5 3.9 335 4.9 45 5.85.5 6.8 RAD 






Pigure 5.24. Measurement 3--Pregquency-Domain Impulse Res- 
ponse (285 mn) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Maesteg. leer i2 2 OIPOLE ANTENNA #12¢28S5NM)> 


MAXIMUM PEAK VALUE. cc ccccccees 428.80 mVU-M 
MIHIMUM PEAK VALUE... cnscceees “896.67 MUM 
RMS VALUE. ccc ecccscccccccsesees %368.85 mU—M 
MEAN UALUE. . ccc cc ccc cccccccece § WS14.14 MUM 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE # 8.5 





9.50 1.00 1.59 2.90 2.50 3.98 3.58 4.90 4.59 3.88 n8 


Figure 5.25. Measurement j--Ramp Optimized (285 mm) 


193 





the 


not 


MEASUREMENT 4--ARBITRARY COMPOSITE TARGET 
1. Description of Setup 
This measur2ment was made to observe the response of 


OS to a completely arbitrary, but axisymmetric target 


encountered by the system previously. The composite 


target was a combination of target 11, a cone, and target 3, 


aeecylinder,. A drawing of the composite is given in Figure 


5% 25. 


oe be Ae 
c ip a 
Pigure 5.26. Measurement 4--Composite Target 


2. Result 


in 


The results are given in Figures 5.26 and 5.27, the 


physical optics shape and ramp responses respectively. The 
corrected diameter is 11% larger than actual, and the 
corrected length ee Shorter than aerial. The 
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uncorrected length 1s 13% greater than actual. The shape is 
a close approximation oof the targets actual shape, although 
somewhat shortened. The cone portion apvears to be about 
onej-nalft of the the total figure, as it was in the actual 
target. 


Bye 


{tQ 


onclusions 





Based on this observation and the nany others taken 
during the implementation of the OS, it can be said that the 
OS can be relyed upon to produce recognizable representative 


Shapes for known or unknown axisymmetric metallic targets. 


195 





TARGET: RUN DATE OLST ANT TGT REMARKS 
Seessiiesi oliver 4 272? COMPOSITE ARBITRARY TARGET 


DIAMETER CINCHES-CORRECTED) 6.5659 ( @.0800) 
WIDTH CINCHES-CORRECTED> 8.90984 ¢ 0.8008) 
WIDTH <INCHES-UNCORRECTED) &% 13.5252 ¢ 9.0880) 
UMBER OF WAVEFORMS AVERACED 2 21 OPTIMIZATION VALUE = 8B 


RAMP RESPONSE -- PHYSICAL OPTICS 


Pigure 5.27. Measurement u--Physical Optics Shape 
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TARGET? RUN DATE DIST ANT TCT REMARKS 
Creaali-Si° le2rn 4 2772? COMPOSITE ARBITRARY TARGET 


MAXIMUM PEAK VALUE. cc cceccens +33.86 mU-M | 
MINIMUM PEAK VALUE. ... -cccvene “335.77 mU-M wy 
RMS VALUE... cc ccccccsevcscccces 4171.95 mmm 
MEAN VALUE. ..ccccncccscccevene =113.28 m=? 
NUMBER OF WRUEFORMS AVERAGED = 21 OPTIMIZATION VALUE = B 





Fae 
ee te | 


Q.50 1.80 1.58 2.88 2.58 3.08 3.59 4.08 4.50 5.98 n$ 


Pigure 5.28. Measurement 4--Ramp Response 
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E. MEASUREMENT 5S--DIELECTRIC CYLINDERS 
1. Description of Setup 

Two runs were sade using different dielectric 
cylinders in place of metallic objects as targets. The 
dielectric cylinders used are similar in configuration to 
cylinder targets 2 and 2 machined from aluminum stock, but 
are made of plastic. For this measurement, the di2lectric 
cylinders will refered to as target 2 and 3, respectively. 

The purpose of the neasurement was to observe the 
response of the dielectric cylinders to axially incident 
electromagnetic fields and to also see how the OS responds 
to a target which was non-metallic. The targets were placed 
fo che standard control location, 1.27 meters from the 
transmitting antenna. The runs assumed the targets were of 
arbitrary, unknown shape for 9S observation. 

2. Results 

Taeget 2 was run first. The results are given in 
Sagquaos 5.238 to 5.31. The physical optics shape in Figure 
5.28 is a very poor representation of the true targets 
shape. The reason for this can be observed in the ramp 
response of Figure 5.29. The S/N is extremely low. A noise 
spike at about 1.6 nsecs has caused the entire resultant 


shape to be foreshortened. The amplitude discrepancy can be 


oS 





attributed to the fact that the amplitude compensator was 
developed for the relatively large valued returns of the 
metallic cylinders. Nevertheless, when the ramp response is 
optimized, the general shape of the target can be seen as in 
meaure 5.30 and 5.31. 

The results using Target 3 are shown in Figures 5.32 
mo 5.35. Because this target is physically larger than 
earcet 2, it has returned a substantially larger signal, 
large enough so that system noise is less of a significant 
factor (although still seen on the waveform). The amplitude 
is again much below that to be expected for a metallic 
eatgst of Simiidiar size. 

3. Conclusions 

Although the information provided by this 
measurement is limited, it appears that dielectric targets 
might provide returns which can be identified by time domain 
techniques. This could be an important area for more 
intensive research as the use of high density polymers and 
melsens in teal targets of interest is of growing importance 
to surveillance systems attempting to locate and identify 


such targets. 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
Seecert-Sh leer 3 722, DIELECTRIC STEERER (3X12 IND 


DIANETER CINCHES-CORRECTED) 0.7662 € 68,8080 
WIOTH CINCHES-CUORRECTED) 2ulo76 C50; 2000) 
WIOTH CINCHES=UNCORRECTED x 3.23813 ¢ 8.8088) 
NUNBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 


a ED ED ED SD SD aD SD ED SD SD SD SD SD SD D> S@® @ ® DP SD SD SD SD SD SD SD SP DP SD SS SS SF DS SD DD DS © ODS S] @ OF DPS FP FOG OPPS @ @ D@ COG SS OHSS Sse 


RAMP RESPONSE -=- PHYSICAL OPTICS 


Figure 5.29. Measurement 2--Physical Optics (Target 2) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 
eabo-t1-St 1.27" 1 27? DIELECTRIC CYLINDER (3X12 IND 


PMN OMUTIOeERR VALUE... ccc cn cscs 11.15 mU-4 
Men PEAK VALUE. sc. cscs cee ° “28.59 mU-M 
Rite VRIEUE . 6 ccc ccc cc ceses oe 417.07 mVU-M 
MEAN URLUE. ccc cccncccrancccnes “13.94 mU--I 
NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 8 











mU-N 
|) Sea a ree ee ae ee 
= ; e 
3 a_i a oe oe a oe ee ae ea 
“5 2, Aa 
-~1or._| 


a Oe a See ee 
isc 1h [Yin ma AU | TT A 
cot LY A, te 
a TE Se a a a 
= | ee ee a 
ee | Ue le 
(i UR Si 2S ee oe ee 


6.50 1.80 1.59 2.88 2.589 3.90 3.59 4.00 4.58 5.98 nS 





Figure 5.30. Measurement 5--Ramp Response (Target 2) 
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TARGET: RUN DATE DIST Alias TGT REMARKS 
Seubert .e ctl ae are SA EOREEE (3X12 IND 
OIANETER C INCHES CORRECTED) 8.°662 8.088 
WIDTH CINCHES-CORRECTED> Ce 1596 ; a. 3008) 
WIDTH <INCHES- ~UNCORRECTED) £3 9.0828 <« 0.0898) 
NUMBER ge _HAVEF ORNS _AVERAGED = 21 OPTINICATION: VALUE = 1 


“RAMP RESPONSE -=- PHYSICAL OPTICS 


——> 


Pigure 5.31. Measurement 5--)ptimized Physical Optics (Tar- 
get 2) 
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TARGET: a aay OIST ANT TGT REMARKS 


seas l 1.274 27? DIELECTRIC CYLINDER c3K12 IND 


Weel PEAK VALUE... ccccccees 48.29 a-i4 
Ripert PEAR VALUE... .. 50.0 ies -27.71 mMU-N 
RNS VURLUE... cc ccc cece cesccens 413.28 mU-M 
PIINMMIUE ss oct cece ec wees eees “10.39 mU-M" 
NUMBER OF HAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE = 1 © 


ee 
| See eee 
ee ee i: 





9.50 1.989 1.59 2.09 2.580 3.08 3.58 4.08 4.58 5.86 nS 


Pagure 5.32. 


Measurement 5--Optimized Raap (Target 2) 
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TARGET: RUN eee oer ANT ial REMARKS 


Wisse l.27 l ee eae ee ones (6X6 IND 
OTA eee cINCHES- CORRECTED) 1.4154 ¢ 
et INCHES-CORRECTED See ; 2 9008) 
INCHES- UNCORRECTED I ¥3 3. - 0089) 
NUMBER OF  WRUEFORMS AVERAGED = 21 SPTIMIZATION VALUE = @ 


RAMP RESPONSE -- PHYSICAL OPTICS 


_ i EPA, 


i iii NS Higa 


Figure 5.33. Measurement 5--Physical Optics (Target 3) 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
Gape-ileGt f.2er" 1 22? DIELECTRIC CYLINDER (6X6 IND 


MAXIMUM PEAK VALUE... .cccccees +94,83 mU-M 
MINIMUM PEAK VALUE... cc cscccece “71.29 mU-{1 
RMS VALUE. .cccrcccscccscvcsace 435.90 mU-M 
MEAN VALUE. .cccccccccccesccces -3.50 mU-M 
NUMBER OF WRUEFORMS AVERAGED = 2i OPTIMIZATION VALUE = 8 


39 








AL 
CeCe 


-30 1.08 1.58 2.890 2.50 3.08 3.589 4.08 4.38 3.80 nS 
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Figure 5.34. Measurement 5--Ramp Response (Target 3) 
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TARGET: RUN DATE OIST ANT TGT REMARKS 
eo G-.1-a1 i.arM 1: 77? =OLELECTRIC CYLINDER (6X6 IND 


MAXIMUM PEAK erase 6 « 6 cisle 6 8.99983 aU 
MINIMUM PEAK UIE cists cers ose 6 -@.9004 nV 
RMS SHEN Y Geclctanovaita.o 6) elese-cxe (es ee eve «© 9.9901 mY 


MEAN oN Ee | a eoeeeeeasv7#edeese?#e 0.0000 mY 
NUMBER OF WAVEFORHS AVERAGED 2 21 OPTIMIZATION VALUE # 9.5 


8.9882 


8.8880 


-8@.8882 
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Figure 5.35. Bea ene S--Optimized Physical Optics (Tar- 
ge 


206 





TARGET: RUN DATE ee ANT TGT REMARKS 
Seti=3 sent 3 7??? DIELECTRIC CYLINDER (6X6 IND 
NUMBER OF WAVEFORMS AUERAGED 2 21 OPTIMIZATION VALUE = 0.5 
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PIPUESE Fe T=—PHAse 


PTT T TTT TIT ty 
Heh A RIO ala DOL bab 
Oe a A a a DC 
° 22 0 eee 


2.9 8.5 1.9 1.5 2.0 2.5°3.8 3.5 4.07425..5.89 5.5 6.9 RAD 










3.88 


8.88 





-3.08 


Figure 5.36. Measurement 5--Optimized Ramp (Target 3) 
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VI. SUMMARY 


A software Operating System has been developed fer the 
Time Domain Radar Laboratory. The System has been desiaqned 
to b= relatively easy for inexperienced operators toc use 
Mec ah Unitial introduction to the OS code and the 


constituent hardware controlled. The OS xtends the 


th 


Sapaotiasy Of the IDRL to the point where accurately 
processed transient time-domain data is presented in a form 
Suitable for operator evaluation. The 0S should serve as a 
DOrimay research tool for the study of electromagnetic 
transients at the Naval Postgraduate School. 

The oss Hh2guLy Operator procsessor interactive. 
Although the code that makes-up the OS is relatively 
Go@plicated, its actual implementation and use is not 
daeimeeeGuit . There are numerous directiv2s issued by the 
System and displayed on the Graphic Systems Crt to the 
operator requesting interactively the input of information 
through the GS keyboard at vital points in the programs. 

The OS is composed of three subprograms. They are: 

laeeeeus for the inputting of Signal information fron 


the ITDRL ground plane 


208 





Gams Maehetai{ical 


} * 


Pee fOr) 6pSeiormang «©6¢the)«6Ssignif 
computations for the 0S 

3. GRAPH for the presentation of the processed 
information in a form suitable for user evaluation 
and analysis 

INPUT receives the desired signal in discrete digitizad 
form from the systems Digital Processing JScilloscope, the 
Peene pat analog to digital converter us¢i in the TDRL. The 
source Signal is a rapid rise time gaussian pulse simulating 
an impulse. The acguired information is averaged a nuaber 
of times as interactively requested by the operator. This 
is done to reduce the white gaussian noise content of the 
acquirsd analog signal. 

MATH performs ice aeiha shag conversion and Signal 
deconvolution to determine the impulse response of the 
target. ihemaiso amplements 42 noise redustion algorithm to 
Optimize the processed data. The amount of optimization is 
under the control of the op2rator and is dependant on the 
degree +9 which loss of small surface features can be 
tolerated and yet obtain 2 good ctepresentation of the 
targets true shape and dimensions. The step and ramp 
responses are also obtained and stored for Later evaluation. 


MATH determines the Physical Optics Snape of the target 
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eRe lFamp C2eSponse wncorrected f5r the effects of 


base2d on 
ig GE. 


p) 


the + 


Owe GS Ody. Of 


induced currents 
and appliss compensation factors 
e rrom 


computes 
historically 


observed 


ee alayal based on 
The compensators are applied directly 


pre-defined targets. 
of 0% Or sbetter 


+o unknown targets resulting in an accuracy 


med termination of the true target dimensions. 
Obs ene 


the processed data to the GS CRT 


GRAPH outputs 
graphed includes the real 


4052 Graphic System. Information 
tine direct, ine ide Me. augmented and backscattered 
waveforms, the time-dicmain ramp crespons?, physical optics 

domain spectral data 


target and the frequency 


Wavertorms as well as the 


Shape of the 
direct and backscattered 


BOL fre 
trequency domain impulse response. Memeans £Cr 2a ote. sng 
the return of progran 


the optimization variable and for 

SomerO. FO LNPUT to aCquire a new Signal is also provided. 
The total result is a highly versatile and powerful 
transient analysis tool. Pi eOseeies  LlLoegnbLlity in ats 
future expansion and/or 


coma LlOworelLoe 
Gonst. vent: 


Pmeesmal Soganization 

Bet racation so aS to allow the contro. of 
hardware and software in applications involving other 
aspects of time-domain target identification. This is 
Met hod 


demonstrated in the complete integration of a Prony's 


ZU 





fyeemed «Solution for the location of a targets poles and 
residues Wine A the requirements Eoz pure Slanocehe 
time-domain analysis procedures. 

The code itself is complests as it stands. However, there 


aemeewceas in Which ge could bea mad= mncre efficient. 


Pecticulazly, a number of graphic subroutines could be 
combined into one routins so as to reduce memory 
reguirements. Also, a method for the display and sélection 


Of menu items using placement of the screen cursor to select 


}4- 


the tem might be developed. This would allow for the 
display of more information graphing processes than is new 
allowed. 

One of the specialized algorithms developed for the OS 
4S a method for expanding certain displaved graphs while 
Gemamneng a2 true plot of relevant information. The result 
is that the GS has essentially infinite resolution and would 


alebowe tor a4 close 199k at particular areas cf the display 


Without disturbing the true graph relationships. Although 


+2 


fameelveteplzed t> Only one SUDrOULIne in GRAPE, it could 


a 


pr 
easily be adapted tt. any other graphic subroutine. 

A second specialized routine is one whereby a true 
one-to-one relationship between coordinate axis is 


Maintained. The unique feature of this algorithm is lis 
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generality in being applicable to anv size date field on any 
Beze Viewoort without ever les2n¢g (mrormaz=On aue ce 
MyiewDOrs clipping". 

The results of field observations using “he 6S have 
proven remarkably faithful to theoretical rasuits. However, 
MeempLesent use of the OS has been principally #irected 
towards verifying that the system was up and working. As 
such, data was collected for axisymmetric targets almost 
exclusively as there was a arjge amount of literary support 
in this research are2a. There still is a need to explore 
System operation with this configuration sf targets. Tn 
Particular, the OS needs to be extended so as to be apnle to 
Bomven ine) Sagnecic rield Integral equation for EM transients 
on the surface of a target. This calculation would then be 
G@oOWe2read With field results to obtain a corrected and 
accurate representation of the target both in shape and 
dimension. 

Hardware has proven very well suited to the needs of the 
2D ab) However, a monostatic antenna more directional chan 
the present omnidicectional radiator presently used, may 
Make the system more generally i2pplicable. Improvements in 
reception by better match of the systems i2ntenna components 


Suowmd allow for the determination of finer feecture on 


2, | 7s 








smaller targets. A more rapid rise tine generator having 
gredtet power outputs than presently supplied by the UHF 
impulse generator would be useful. a naa, ays fulay 
computerized analog digital processor with a microprocessor 
controlled sweep rate and distance windowing function would 
eagommucy tO the automation and, consequerntially, consistancy 
Se cosults in the TORL. 

pid le project is dynamic and holas promise £OE 
potentially practical applications of applied tresearch. 
Ther= rémain many areas t9 explore by wsing the 05 
lmaginatively. The system is now ready. The true need is 
for the refinement of procedures in using the OS and 
applying it to research areas having merit both instructive 


ama. practical. 
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AS stated in Chapter IV, three peripheral ROM packs nave 
been added to the Graphic System to tailor it to the 
specialized needs of the MTDRL. iiey Jabe- che Sitonal 
Processing ROM Pack No. tT, the Signal Processing ROM Pack 
Nee 2 (261) and the editor. 

ives gaat Processing ROM Packs provids fiftesn 
Mathmatical operations that are essential to advanced signal 
precessing and data analysis. Suge processing is 
performed on waveforms or data stored in one-dimensionsal 
azrays. These capabilities do no alter the operation o= the 
GS or occupy any of its avaliable RAM space. A summary of 
each ROM command follows. [Ref. 24 25} 

A. Signal Processing ROM Pack No. 1 
1 MAX -- finds the maximum of an array. 
Coie eee tnds the Minigum of an array. 
3. CROSS -- finis the location of a specified crossing 


level within an array 
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Beez s- PerrormS a2 two-point dirterantiation of an 


an 


ri) 


Dera PeSaEOrmS Aue cnree=point ditferentiation o 
array. 
Mie == Ineeqrates an array. 
BEoey—— dssetays a graph of. an array in raw <#orn 
(without axis). 

ama Processing ROM Pack No. 2 
Patees=) COmputes the fast Fourisr transform of an 
array, placing the results in the san2 array. 
IFT -- computes the inverse Fourier transform of an 
array, placing the resultS in the same array. 
CONVL -- convolves two input Eravsy Veclkacing .) tne 
Rosters an 4 third array. 
CORR -- correlates two input arrays, placing the 
results ina third array. 
Bowe =—- “COnverts sgt ee e@el eels So iaig ah data from 
rectangular form (reals and imaginaries) +9 vpolar 
zorm (magnitude and phase). 
Meek Multa plies an array by a cosine window of 
program-selectable tapering weighs. When tapering 


is selected to be 50%, it provides a Hanning window. 


(AS: 











femeiigbay == SOLtTS an array Of interleaved FFT data into 
mao, 2Ecays, Deoweenctainard Leal and One containing 
imaginary components. 

8. INLEAV -- interleaves the real and imaginary data 
from two input arrays into a third array whose format 
is acceptable to the IFI command. 


irmware 


rth 


Mromer Deon 2S a ROM device that contains 
routines that allow ASCII magnetic files to be altered or 


edited off-line. [Ref. 26] 
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Ae GRAPHING VARIABLES 


The following named variables, numeric and constant, are 
used in more than one subprogran. A brief description of 
mets DpUrpOSS and the file in which thsy ara initiated is 
given below. NOG-ows pect fic PLOgram/7subtoveine for local 


numeric variables, constants and arrays. 


en ee es eg 
| | 
{VARIABLE | PURPOSE BLS 
NAME | | rurrratzp| 
a nr 


| | STRING CONSTANTS | | 


—-----—-+—--— 
2 | 








| 

| BS | DAT Seer RUN 
[+ ———-____-— 4-4 
| 2s | TARGET DISTANCE ( { 

| | 
| GS ANTENNA NUMBER Z i 
HS | TARGET NUMBER | 

| | 

MS VGRPLCALYSCALe aeTOR = mVOLES | 
| 

PS | RUN NUMBER 


| 
| | 
a Lrt—“—s—SC—s—s—sLC 


Al i 





Ee eens 
| vaRraBLe | PURE OS & Ss 

| 

{ 


| 
NAME | 





| NUMERIC CONSTANTS AND VARIABLES 











| 
| 
| | | | 
| 36 OPTIMIZAIION CIOAPENSATION { 3 { 
| | | | 
{ 37 | OPTIMIZATION COMPENSATION | 3 | 
— | | 
B8 | OPTIMIZATION COMPENSATION { 3 | 
| | | | 
| ECS { TARGET LENGTH | 2 | 
[nn 
C6 TARGELD RADIUS | 2 | 
i 
| Co | ANTENNA PARAMETER AVERAGING VAL { Z i 
| | | | 
L { OPTIMIZATION FACTOR { 3,4 | 
| | | 
| ia | FLAG 3 7 
) 
——— 
j NZ | SISNALS AVERAGED/SWAVSFORM { ie 4 
| | | | 
| P7 { FLAG 


OE 


RUN NUMBER 


TU 
(o@) 
tO 


¢ 
2 
| 








ew _——— ee ee eee ee eee ss ee 























cr— 
[VARIABLE | CURE Uss | Pi Le | 
| NAME | eae 
a 4 
| oe { OPU Tut ok beON COMPENSATION | | | 
a 
L a9 | ENDPOINT DISTANCES OPTICS i 3 | 
| | 
S1 ViGemitess ocALS FACTOR | 2 | 
a ne a pe ee ee 

| 
| 
i aaa FLAG | _ 
 ___ 
35 fig a LF oes EL | { 

| ANTENNA PARAMLSTSE ENGTH 
: se 
| | ! 
S6 ANTENNA PARAMETER-WIDTH 3 | 
| | 
i ea DPO SWEEP | Ze 
. | <i 
| i | ———— oe 
a | ies GaGCONSTANTS ARRAY | ape | 
pe 
| Z1 SAMPLING RATE DPO | 2 | 
cee a 


Zag 








A. 


GENERAL DESCRIPTION 


le 


De 


EeGaetOns rile 2 

Memory space required {less nass SeOLEATe). 
24552-bytes 

Maximum no. of steps: 272 (212 less teamarks) 
Description: INPUT provides the means to initialize 
and drive the ITDRL Operating System. It allows input 
of all waveforms, averages acquisitions and stores 
data. ees tows (fOr input of antenna target 
parameters and initializes most jJraphing parameters. 


Inputs are from the k2yboard of the GS, mass storage, 


and File 4, Output 1S to mass storage and File 3. 
This program is Dagny Sperator-processor 
interactive. Numerous messages and directives guide 


the user. 
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Be. LOCAL PROGRAM VARIABLES AND CONSTANTS 


1. Significant 
nr 5 SSIS a, SSS SS == 9 (EN tA 
lyDs | NAME |USE | REMARKS | 
] | 
a 
| ARRAY {3 | DIWAY | ($12) --TEMPORARY Berne Oi 7 
| | | | BIR cer enwavercons 
iE a eee 
| ARRAY re |AUGHAY Me ee On Any ARRAY £ OR | 
AUGMENTED WAVEFORM 

SS a ae cae 

| ARRAY X0 lACAYS | (512) --TEMPORARY ARRAY FOR | 
| | | NWAY SIGNAL ACQUISITION | 
| loraay | AND AVERAGING 
| | ;AUGKAY | | 
4 | {\SCATW AV | | 
ES Se ee | 

| CONSTANTS | 20 |ALL | MENU LTEM | 
| | ee Sara 


See ometNG CONSTANTS: As, 85, CS, 05, JS, L3, KS 


C. NUMERIC VARIABLES AND CONSTANT 


U2 
©) 
=a 
= 
O 
tO 
O 
Ud 
Q 
4= 
~ 


Ove NaN P9ees0, T, x 
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C. MEMORY REQUIREMENTS 

















Eo ocelleernal Daza Scorage 

Ss LS dUhUmUmUmC”t~—<“‘“O™SOSOSOCOCQNNO®”™”;:;‘SQ 
DES Gaeet LON ok) Gallntyee’ { 

a | | 

| pacH ACQUIRED WAVEFORM TO BE AVERAGED (4114 { 
| | | 
{AUGMENTED WAVEFORM ARRAY (4714 { 

| 

ae ee ) 
| prRect WAVEFORM ARRAY ata { 
oe Ee ae 4 
| yaRroUS NOM ERICSSTRING VARIABLES AND CONSTANTS [5575 i 








| 
SOTNE tT ot ee 


--- 





—— 


MAXIMUM MEMORY REQUIREMENTS (INTERNAL) FOR INPUT 442869 




















2. External Mass Storage Reguirements (Data Taps) 
are ae ee a 
| DESCRIPTION | BYTES| PILZ | 
ae | i 
| prRECT WAVEFORM ARRAY eer | 50 | 
| INCIDENT WAVEFORM ARRAY lutts | 51 
a ane er: ean 
| AUGMENTED WAVEFORM ARRAY eer | ae | 
EEE Ee a 
| SACKSCATTERED WAVEFORM ARRAY luqtg | 53 | 
ao a a ower 
Le TOTAL lan | 
i EEE pe ee 





SUBROUTINES 


D. 




















=) > fe] m | 
mM © OQ fry = = 4 
tS 4 an > > nN uy «tt E> «cd 
tH E44 QO YN fc (4 fr] et wt E+ gq €4 
| [| ti Ei c= oo << eG = +t 
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C) fx ws A ey Pos me | = uw 
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mm wU? Oy at oy) ee S| G4 = E+ fr, fx 
nm WwW fy E4 tr Q) » ) ~ I hmQy Oo % 
a cH MO HH DD OM H WW fu Ww S 
oO Bb} @& H OA mh es Ph mm hr oO 04 . 4 
Hi < tty ~ iz, | MoH ta > ah W 
fA - ° mm WM KF HM FI H Mier Ei fits 
D4 Ys oH | mM ww ra 2 | =) Hi 
4 Lad bf |) Coamned cls Mee xm Wh mr UY) ct oM wy Aa > 
03 tea fx | HH et ee | Hi t{ fr | & at 
t) Fa tn [WES 364) Le Fa 0G oo ee 8 or He | at “Sto ore 
un) HH oO Mm f2 W Ke © Ft! Oo tH OO Or DO 65 
a a Mm O OO & PD in a ZZ = = 
Gq Ma Oy £4 xq mf =r eM Ww Me Ww HW Hifi «ft 
md .O 
-—~ = = ee SP ap Se eee SF a cas “OEP ee 
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Zo 
| 
6 
6 





9 





‘s 
| 


—————— 
JLrg 
+----—+ 
aro 0 
ins 
—__}___- 
| >000- 
‘an 
i— 
{ 3000 
[3170 
———}-—--—-| —---+---}L —---- 
{19 
| 


| 3200- 
13380 
ae 


t 
| Se Ir Oss ] 
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}—— - + —- + —--- —} -— -~ = --- 4, - ——+ 
CLASS [was [LINE | STEPS 3 | DESCRIPTION | 
| | 

pan nf 
ISPECIAL{ANTLOC {| 4000- 3 {| 2 {INPUTS ANTENNA PARA. | 
L | | 4060 | | | FROM MASS STORAGE : 
——-f-—-- =f -— ---} + -- --- 4 

| spzcraLlt TGILOC ae \ 7 > PaAnGel Dew eENSLONS | 
| peek ieee eons FRO INA SS | 

| | iSTORAGE | 

| | | 
SPECIAL|NEWPAR [4200- 4 1 | CNITIALIZES ANTIN AND 


a 
E 








| | 
| | 4 | 
| | | 4330 TGTIN FOR STORAGE OF | 
{ { | | NEW ANTENNA AND 
i i ITARGEI PARAMETERS | 
soos (een eran al Pa Pee 
{|SPECIAL{ANTIIN j4400- {21 | 7 {INPUTS NEW ANTENNA | 
| | | | = a - | 
[_ | {42600 i a aa EO “Aso 7 
| | ! 7 | . | 
- | ero [it | ° -_ # NEW TARGET | 
4800 i i {PARA. TO MASS STORAGE] 
eee ee 
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GENERAL DESCRIPTION 


BocaciOn: File 3 

Memory space reguired (less mass SeOrage),: 
14544-bytes 

Maxinun no. of steps: 206 (169 less remarks) 
Description: MATH perrorms oie Signa f2ecane 
mathematical computations for the operating system. 
These are Ease DOUELeE and inverse HOU et 
transformations, determination of inpulse, step, and 
ramp responses, noisé optimization and compensation, 
and antenna parameter determination. Inputs are from 
File 2, File 4 and mass storage. Output is to File 4 
and mass storage. ess sees Cared 1 1s competely 


automated with no Operator-processor interactions. 
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A. 


APPENDIX 2 


GENERAL DESCRIPTION 


Ms 


Ze 


mocataon. File & 


Memory space a-qulLced (less nass storage): 


4&u712-bytes 


Maximum no. of steps: 584 (473 less remarks) 


Description: GRAPH provides the means to visually 
display the results of signal 2eCUlSse2om 2nd 
processing. Subroutines are divided into three 
classes: 

a) GENERAL~-those subroutines that are2@ applicable to 


D) 


most of the other routines in GRAPH. They 


initialize the progran, POENace Graphics, citi 


(D 


v 
lable, and overlay a grid. 

INVERSE SCATTERING--those routines applicabl2 to 
the inverse scattering solution of time domain 
electromagnetic transients. They are the 


secondary drivers for the develooment of graphics. 
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and mass storage. Outputs are to the GS CRI, mass storaye, 
meprce 2, FLis 3, File 6, and File 8. Tavs preogaam 15 highly 
Operator-processor interactive. Numerous messages and 
directives guide «he user. Signicant graphs are provided 


for operator evaluation of processing results. 
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In this Appendix is presented a typical run in which a 
Sphere and a con2 are llluminated by The ine tdene 
zlectomagnetic field. The messages oresentsd and the 
results displayed are those as would be encountered during a 
normal run of the systen. 

iemOsom as ElrSt Weed £5 a@gquire and display a normal 
Signal (the sphere). Pais Segnal- | tS “then OPtTimized to 
remove noise. The results of the effected waveforms are 
then represented. Finally, the system 1s returned to INPUT 
and a new target, a1 cone is acquired. The results for the 
cone are presented. Note that the impulse response is 
optinized by a factor of .5 while all other graphs for the 


con2 are shown unoptimized. 
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